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ABSTRACT 



An internal airblast from a missile explosion produces tearing, 
holing, and dishing of a ships' structural plating. 

In order to analyze the residual strength of a ship after an 
internal explosion, a typical ship compartment was considered to be 
composed of flat square plates subjected to uniaxial compression, and 
the damage effects were examined independently. The holing was 
considered to be a centrally located circle. Plastic limit load 
analysis of a flat plate indicates that ultimate strength decreases in 
proportion to increasing hole size. Elastic-plastic analysis indicates 
that ultimate strength is little affected by hole size until the holing 
reaches half the plate width, and is essentially equal to the ultimate 
load described by the effective width concept. An extension of the 
concept of small imperfections to dished plates suggests that the 
ultimate strength is reduced by roughly 10 percent from that of a flat 
plate. The extensional collapse mode analysis describes folding and 
strength beyond peak load. Fully plastic crack propagation was 
examined, but its limiting effect can not be precisely described. 

A small program of experimentation was developed to gain 
qualitative insights into the interaction and cumulative effect of the 
various damage modes. Seven mild steel boxes were constructed of four 
square side plates with length- to- thickness ratios of 90 and two 
stiffened end plates to model different damage effects. The 
experimental results indicate that the plastic limit load does not 
define the peak load of slender plates. Increased hole size appears to 
decrease ultimate strength gradually . Dishing on the order of 5 to 15 
plate thicknesses reduces ultimate strength by roughly 15 to 3_5 percent . 
The mean crushing load allows damaged structural plating to continue to 
support loading on the order of 5 to 10 percent of the tensile yield 
load. Although crack growth was not a factor during testing, it could 
be important in shipboard structural plating which is relatively thinner 
and harder than the tested specimens. 
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NOTATION 



Most notation will be specifically defined where it is first 
introduced. The following table will serve as a reference guide. Where 
more than one meaning has been assigned to a symbol, the correct use 
will be obvious from the text or explicitly defined. 

A Cross-sectional area, crack area 

a Length of plate 

b Width of plate 

b ,b Effective width of unperforated plate, corrected 

e © c 

effective width for holed plate 

C Coefficient of deflection function 

r 

Diameter of hole 

Width of tensile region in unperforated plate 
Plate flexural rigidity - Et 3 / [12(1 - v 2 )] 

Material constant 
Young ' s modulus 
Airy stress function 
Half-width of the folding zone 
Kinetic energy 
Plate thickness 
Impulse 

Plate buckling coefficient 
Wavelength parameter 
Stress resultants (membrane) 

Material constant 
Unit stress resultants 
Load 

Critical load of unperforated plate 
P^ Compressive load removed by hole 

P Mean crushing load 

m 

P^ Tensile yield load (squash load) 

Q,q Generalized stress, strain 

R Fracture toughness (essential work of fracture) 

r Coordinate, radius of circular plate 6c tube 



c 

c 



E 

F 

H 

h 

I 

K 

m 

N , N ,N 

x y xy 

n 
n 
P 
P 



c r i t 
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\ 



u » u , > u 
b m 

U,V,W 

V 

V 

w,w p ,w c 



w , w 

c c f 

x,y,z 

Y 



a 



P 



e 



e , e , e 

xx y y x y 




e 



r 



A, A 

r 

Ap 

v 

P 

a 

a ,o ,a 

XX y y X y 

a 

ci it 
O 

o 

G 

Y 

O 




Time 

Strain energy (bending, membrane) 

Displacements 

Volume 

Velocity 

Work (plastic, crack propagation) 

Central deflection, final central deflection 

Structural coordinates 

Flow strength 

Johnson's damage number 

Slenderness ratio 

Extensional strain 

Plane strain components 

Average hoop strain 

Effective strain 

Total plastic work per unit volume 

Edge displacement parameter, deflection modeshape 

Elastic strain energy (residual) 

Reference slenderness for plate buckling 
Poisson's ratio 
Material density 
Normal stress 
Plane stress components 
Critical stress 
Flow stress (damage stress) 

Yield strength 
Effective stress 

Modified damage number, circular & quadrangular plates 
Distance parameter 

= d\ )/9x A + 2 d\ )/dx Z dy Z + d\ )/dy‘ t 
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1. INTRODUCTION 



1 . 1 The Antiship Missile Threat 

Modern naval warships must be built to survive and operate in a 
multiple hostile threat environment. In addition to facing heavy seas 
and storms, a warship may be subjected to attack by a variety of weapons 
such as torpedoes, mines, gunfire, bombs, and missiles (Figure 1.1). 
Some of the more common and lethal weapons are antiship missiles 
equipped with semi-armor piercing warheads with small, high explosive 
(HE) payloads. Such antiship missiles first penetrate a ship's outer 
structure and then, by a delayed- fuze mechanism, detonate within the 
confines of the superstructure or hull to cause damage to the structure 
and internal systems [1], The purpose of this paper is to analyze the 
structural ability of a destroyer-sized warship to withstand the effect 
of explosive loading associated with a missile attack. 

The various antiship missiles produced by different nations all 
have similar components. A typical antiship missile has a propulsion 
section (multi-stage solid fuel and/or a turbojet engine plus a booster 
for surface/submarine launch), fuel tank, guidance system, and a warhead 
enclosed in an aerodynamic cylinder. Small antiship missiles are widely 
produced and exported throughout the world. Their sizes range from the 
now infamous 17 foot French EXOCET with a 200 pound warhead down to the 
10 foot Norwegian PENGUIN with a 264 pound warhead. Most antiship 
missiles attack their targets with a sea skimming trajectory which make 
them difficult to detect and track. 

Modern antiship missiles' blast and fragmentation warheads are 
little more sophisticated than the free-fall bombs used against ships 
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since the start of World War Two. 



The effect of aerial bombs and 



kamikaze aircraft attacks on destroyers during World War Two illustrates 
the type and extent of damage produced by such weapons. A direct hit of 
an aerial bomb causes damage to the above -water and underwater sections 
of the hull (sometimes with a breakoff of an extremity) , fragmentation 
damage to exposed stations and superstructures, flooding of compartments 
producing a change in draft, and the outbreak of fires. The main 
propulsion spaces, electrical equipment, and the ship guns could be put 
out of action as a result of a high explosive's blast and fragmentation 
damage. Aerial bomb explosions put destroyers out of action basically 
because of damage to guns and technical equipment (about 50 percent of 
the cases, of which almost 25 percent were due to shock vibrations) and 
because of loss of watertight integrity (35 percent). Fires and 
explosions and destruction of the overall strength of the hull of the 
ship had less significance. The explosive action of kamikaze aircraft 
most often caused destroyers to lose mission capability rather than to 
sink [2]. Of 30 U.S. destroyers sunk during World War Two due to 
gunfire, bomb or kamikaze attack, only six resulted from loss of 
structural integrity of the hull (Figures 1.2, 1.3) [3]. 

Significant ship structural damage due to weapons effects has been 
scant since World War Two and has generally not been as well documented 
[3]. During the recent Falklands Conflict, however, the British Royal 
Navy faced a multiple hostile threat environment. HMS Ardent , a 3250 
ton frigate suffered four aerial bomb hits before it was "mission 
killed" and required three more hits before it sank. HMS Antelope , 
sister ship of HMS Ardent, sank after one 500 lb. bomb exploded and 
started an uncontrollable fire that caused several secondary explosions. 
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HMS Sheffield was hit by an EXOCET antiship missile which failed to 
explode, but which penetrated and ignited the ship's main fuel supply 
tank with its still burning rocket fuel. The ensuing fire and smoke was 
uncontrollable. The 15,000 ton roll-on/roll-off container ship Atlantic 
Conveyor was also lost due to a single EXOCET missile hit which produced 
uncontrollable fires [4] . 

1 . 2 Modern U , S . Navy Warship Design 

United States warship structures have always been implicitly 
designed with survivability in mind. Incorporating lessons learned from 
World War Two into such design practice as crack arrestors, requiring 
gradual transitions, using continuous welds, providing hull girders with 
a 1.0 ton per square inch reserve primary strength and limiting design 
primary stress throughout the structure, has produced more survivable 
warship designs. David W. Taylor Research and Development Center has 
been developing material fracture control technology. Strict 
fabrication requirements minimize the deleterious effects of poor 
structural detailing, failed welds, stress concentrations, etc. in high 
intensity dynamic pressure environments and also contribute to the 
survivability of the warship [3]. Current techniques for predicting 
weapon effects on structures are, however, limited. To simplify the 
determination of residual strength after an attack, the Naval Sea 
Systems Command considers all damaged structure to be "ineffective" [3]. 

Small, unarmored, warships such as destroyers are likely to suffer 
local structural failure from the initial explosion and be faced with 
spreading fire causing further damage and secondary explosions. The 
large amount of flammable material on modern warships make them more 
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susceptible to fire damage than World War Two warships . One of the 

lessons learned from the Falklands Conflict, and further emphasized in 

the Persian Gulf, was that aluminum structures failed because of the 

intense heat produced by fires. Although steel structures may be 

o 

weakened or distorted by shipboard fires on the order of 1800 F to 
2000°F, aluminum structures will rapidly sag and melt at such 
temperatures [5]. Since World War Two, small warships were built with 
aluminum superstructures to reduce topside weight and associated costs. 
The U. S. Navy's new DDG-51 Arleigh Burke Class destroyer incorporates 
improved design practices, all steel construction, and damage control 
features such as shock mounted and blast hardened structures and 
equipment along with redundant firefighting capabilities. Such warships 
can undoubtably survive the secondary effects of missile damage better 
than the the destroyers of World War Two or the British warships in the 
Falklands Conflict. 

Contemporary tactical doctrine calls for coordinated saturation 
attacks by antiship missiles to overcome a warship's active defenses and 
ensure several missile hits [4] . Even with severe battle damage, a 
hardened structure provides a warship with a "keep afloat" capability. 
To improve this passive defense capability against the growing antiship 
missile threat it is important to understand how a warship's structural 
plating is affected by the damage produced by an internal explosion. 
The capacity of a warship's structure to absorb and contain missile hits 
is vital to its ability to survive and operate in a hostile environment. 

1 . 3 Local Damage Effects and Method of Analysis 

An internal airblast from a missile explosion imparts reverberating 
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high- intensity pressure pulses against ship structural plating. The HE 

o 

detonation produces gasses with a temperature on the order of 5500 F and 
pressure reaching 750,000 psi. The total impulse is a function of 
missile momentum and explosive energy. After passage and decay of the 
initial shock wave, confinement of the explosion gases leads to a 
pressure loading with a typical duration on the order of several 
seconds [ 6 ] . 

Structural damage is caused by a combination of (1) holing of 
internal structural plating, (2) blast overpressure within confined 
spaces, (3) fragments from warhead break-up, (4) fragmentation/spalling 
of interior bulkheads and decks, and (5) fire [1]. These modes of 
failure interact and may cause structural collapse of a local or global 
character. The effect of fragmentation on structural integrity is 
generally insignificant and is not evaluated. The subsequent effect of 
heat and fires produced by missile damage on structural integrity has 
been analyzed by Hiddemen [7]. The structural effects of supporting 
stiffeners (girders, beams, etc.) have been examined by Faulkner [8] and 
others. This paper is specifically concerned with the residual strength 
in a ship's structural plating which has suffered an internal airblast. 

The dynamic interaction of damage effects within a complex ship 
structure is not well understood. In order to analyze this situation, a 
typical warship compartment structure is considered to be composed of 
flat square plates subjected to quasi-static uniaxial loading, and the 
damage effects are examined independently. Structural damage consists 
of some combination of tearing, holing, and deflection. 

The extent of holing caused by the missile explosion is based upon 
U.S. Navy experiences [1]. Holing is caused by the high 
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velocity gradients across the plate surface (due to localized pressure 
loading) [6]. Roberts and Azizian [9] have shown analytically and 
experimentally that hole shape has little effect on structural strength. 
As a reasonable simplification, therefore, the holing is considered to 
be a centrally located circle. 

As an initial approximation of residual strength after damage, a 
plastic limit load analysis and elastic-plastic buckling analysis are 
developed for a flat plate with a circular hole. The dishing of a flat 
plate under impulsive loading is examined to determine the relationship 
between lateral deflection and blast overpressure. The residual 
strength of a dished plate is examined by an extension of the concept of 
initial imperfections and by the recently developed extensional collapse 
mode analysis. The limiting effect of mixed mode crack propagation upon 
residual strength is then examined. This progressive analytic approach 
is presented in Figure 1.4. 

A modest program of experimentation is conducted to determine the 
relationship between deflection, holing, crack propagation and residual 
strength. An analysis is then made of the residual strength in a modern 
U.S. Navy destroyer (DDG-51 Arleigh Burke Class) hit by one antiship 
missile amidships . 
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Figure 1.1: 



Categorization of 



naval antiship weapons 



[3] 
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DD-518 
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26 DEC. 1943 
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NAVY DEPT. 


BUREAU OF SHIPS 

- - 



Figure 1.2: Diagram of battle damage to USS Brownson (DD-518) which 

sank due to loss of structural integrity [3] 
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SHIP 
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SIMS 
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\ C*. . 


2100 TONS 
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WITH BOMBS. 


WITHIN 
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Figure 1.3: Summary of battle damage to several World War Two U. S. 

destroyers which sank due to loss of structural integrity [3] 




Figure 1.4: Analysis of battle damage, (a) Real damage (b) Uniaxial 
compressive loading of a flat plate with a central circular hole (c) 
Uniaxial compressive loading of a dished plate (d) Uniaxial compressive 
loading of a dished plate with central circular hole and initial cracks 
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2. PLASTIC LIMIT LOAD OF A FLAT PLATE 



2 . 1 Introduction 

Shipboard structural members such as bulkheads and decks can be 
modeled as thin plane sheets of homogeneous material. Such a structural 
member is often acted upon by forces in its plane which are much larger 
than the normal and shearing forces acting on its plane. If the sheet 
lies in the x-y plane, the non- vanishing stress and strain components 
are o , o , o and e , e , e respectively. The loads which cause 

xx yy xy xx yy xy 

plastic collapse (yield failure) in a structure can be analytically 
described through the use of lower and upper bound theorems. 

Structural damage consists of some combination of tearing, holing, 
and deflection. The compliance of a ship generally imposes a uniform 
displacement (which some authors describe as constant strain) boundary 
condition on the loaded edges. A rough upper bound limit load analysis 
of the effects of fully plastic fracture in a dished plate with a 
central circular hole is presented in Chapter 6. As an initial analysis 
of the effects of damage, however, a plastic limit load analysis of a 
flat plate with a central circular hole subjected to uniaxial 
compressive loading is presented here. 

2 . 2 Physical Basis of Plasticity 

Plastic deformation occurs when a structure experiences loads which 
produce stresses above the yield strength. The stress - strain curve for 
mild steel, shown in Figure 2.1, illustrates actual versus perfectly 
plastic deformation. Microscopically, plastic deformation consists 
primarily of shear relative to crystallographic slip-planes, with no 
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volume change. Three consequences of the micromechanics carry through 
to macromechanics: (1) incompressibility, (2) current stress is related 
to the increment of strain, not to the total strain as is the case with 
elasticity, and (3) the associated flow rule. A continuum theory of 
plasticity can only be applied to relatively large regions as compared 
with the mean free path of dislocations and the grain size (0.001 to 
0.1mm), or perhaps to coherent slip bands. 

There is a relation between the yield locus, which determines the 
state of stress at which plastic flow can occur, and the equations 
relating stress and strain increment for the subsequent plastic flow. 
The associated flow law establishes a one-to-one correspondence between 
the state of stress and strain rate if the yield locus is smooth with a 
finite curvature at each point. It has been shown to hold not only at 
the microscopic level but for forces and moments on structures as a 
whole. Once the yield strength has been reached, the current state of 
stress determines the ratios between the strain increments but does not 
define the strain. A hole in a tensile structure under plastic loading 
will, for non-hardening plasticity, localize the strain in a narrow band 
and that band will extend all the way to the nearest boundary of the 
part [ 10] . 

2 . 3 Limit Load Analysis 

As a structure experiences general plastic deformation, its 
load-deflection curve becomes nearly horizontal. The limit load is the 
horizontal limit approached by the load- deformation curve. Although it 
can be very difficult to analytically determine the limit load of a 
structure, it is comparatively simple to determine upper and lower 
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bounds to the limit load. 

According to the lower bound theorem, any load which is in 

equilibrium with a stress field such that the yield condition is not 

violated and the equilibrium equations are satisfied is a lower bound on 

the collapse load. The equations of equilibrium for plane stress are: 

do do do do 

” + — - 0 — ^ = 0 . ( 2 . 1 ) 

dx dy dx dy 

There are several ways to express the yield condition. Hodge [11] used 
the Tresca yield criterion which can be expressed as: 



max 



O +0 
x x y y 



yT-“) 



+ o 



xy 



rO - O 


, 2 


xx yy 


2 




+ O 


1 2 J 


I xy 



< O 



(2.2a) 



and Gordon and McCrum [12} performed a limit load analysis using the 
Mises yield criterion which can be expressed as: 



/ \ \ \° - o ] + o Z + o 2 1 + 3o 

y / L l xx yyj ** yy J 



< o 

xy y 



(2.2b) 



where o^ is the material yield strength. In constructing lower bounds, 
it is frequently useful to consider discontinuous stress fields across a 
plane with normal n and in-plane coordinates s,t. Then the components 
acting across the surface o f o ,o must be continuous for 

n n ns n t 

equilibrium, while the in-plane components ° tt '° > a t need not 

An upper bound on the limit load can be found by considering any 
incompressible velocity field which satisfies all velocity boundary 
conditions. It is frequently useful to consider discontinuous velocity 
fields when evaluating upper bounds, which is permissible as long as the 
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work associated with such a discontinuity is accounted for. The yield 
locus is a contour of constant value of the flow function which 



characterizes the yield condition. Analysis of the yield hexagon 
reveals that the internal energy is the product of the yield stress by 
the absolutely largest of the three principal strain rates. Any load 
which produces an external energy rate in excess of this value is an 
upper bound on the limit load. For the case of uniform traction (which 
some authors describe as constant stress) , the upper bound to the limit 
load P + can be expressed as: 

+ 



<5W 



(2.3a) 



(l/A) Jsu p dA 

where A is the gross cross-sectional area, 6W P is the increment of 
plastic work, and 8 u p is the increment of deflection along the line of 
loading. For the case of uniform displacement, the upper bound to the 
limit load can be simplified as: 



p + = 



sw 

8 u E 



(2.3b) 



2 • 4 Plastic Limit Load of a Holed Plate 

The limit load of a flat square plate with a central circular hole, 
subjected to uniaxial compressive loading, can be determined within 
reasonable upper and lower bounds. Hodge [11] analyzed the tensile 
limit load of a flat square plate with a central circular hole, simply 
supported edges, and uniform traction boundary conditions. Figures 
2.2a,b,c demonstrate that similar but orthogonal compressive or tensile 
loads will produce the same rigid-plastic failure. This can also be 
demonstrated mathematically by changing the signs of loading and 
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associated deflection in Eqs . 2.2 and 2.3. The rate of external work 

required to produce similar deflections is the same for both compressive 
and tensile loads. 

Hodge constructed a kinematically admissible velocity field by 
considering a sliding of the slab out of its plane along a plane through 
the axis of symmetry. From the Tresca yield criterion, an upper bound 
P + on the limit load of a flat plate of width b and with a central 
circular hole of diameter c is: 

p + - p y a - f ) (2.4) 

where P is the tensile yield load defined as: 

Y J 

P - a A . (2.5) 

Y Y 

0 

For reasonably large values of — , however, a better upper bound value 
must be based upon the bending model, utilizing the concept of a plastic 

0 

hinge. If g is almost 1, then the holed plate could easily collapse by 

means of yield hinges at its four narrowest points. The yield-hinge 

concept can be used as a guide to obtain a reasonable velocity field. 

1 c 

The upper bound limit load for ■=■ < r- < 1 thus becomes: 

5 b 

p* - / 2 • 4 E + 3 <5> 2 • E ' < 2 ' 6 > 

A lower bound for this problem was obtained by considering the 

discontinuous stress field shown in Figure 2.3. The lower bound limit 

load for 0 < £ < 0.443 is: 
b 

p- . b (L: c/b)(l - P ( 2 . 7 ) 

c 

where f is a distance parameter chosen to reflect a discontinuous stress 
field for maximuum P . For larger size holes, with 0.443 < r- < 1 , the 
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lower bound limit load is: 




2 



( 2 . 8 ) 



The upper and lower bounds are combined in Figure 2.4 to describe 
the plate failure loading for different size holes with uniform 
traction. The upper and lower bounds derived by Gaydon and McCrum [12] 
using the Mises criterion, presented in Figure 2.5, are very similar to 
those derived by Hodge. Gaydon [13] also performed a biaxial limit load 
analysis of this case with roughly similar results. 

For uniform displacement of the loaded edges, which more closely 
approximates shipboard structural boundary conditions, limit load 
analysis with Tresca criterion turns out to provide an exact solution. 
The upper bound displacement field is described by two 45 degree 
sections across the mid-plane, as shown in Figure 2.6. The lower bound 
is described by the constant stress regions outside the central hole, as 
also shown in Figure 2.6. The lower and upper bounds coincide to define 
the limit load: 



This simple linear relationship, which corresponds to the upper bound 
for small holes, is depicted on Figure 2.4 along with the uniform 
traction upper and lower bounds. Figure 2.4 suggests that shipboard 
compliance, which restricts non-uniform displacement of structural 
plating, improves the strength characteristics of the plating. 

A separate method of determining collapse load consists of simply 
evaluating the stress experienced by the net plate structure at the 
mid-plane and comparing it to the yield strength for the material. 
Stress concentration effects, in general, apply only to elastic behavior 



P - P 



Y 




(2.9) 
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and are not considered in this approach. In the case of a flat plate 
with a circular hole, the stress experienced in the remaining structure 
is increased and the uniaxial load required for yield is correspondingly 
decreased as: 

1 - gj - (2.10) 

This simple "engineering" approximation is identical to the uniform 
displacement limit load. 

The limit load theorems are based on current geometry. Structural 
deformations may significantly influence the boundary conditions and 
equations of equilibrium, producing discrepancies between the actual 
limit load and that based on the original shape. Ductile structural 
materials may not behave in a perfectly plastic manner. Thin plates . in 
particular, can experience considerable deformation in the elastic range 
which distorts the structure from the assumed plane shape prior to 
plastic collapse [14]. Despite these limitations, if elastic -plastic 
buckling has not occurred, plastic limit load analysis provides a 
reasonable first look at the effects of holing. 
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Figure 2.1: Graph of typical stress versus strain for mild steel [11’ 
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Figure 2.2: Square plate with a circular hole, (a) Tensile load (b) 

Compressive load (c) Uniform traction hinge type mode of collapse [11] 
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Figure 2.3: Discontinuous stress field of square plate with a central 

circular hole (uniform traction) [11] 




Figure 2.4: Graph of bounds for a square plate with a central circular 

hole using Tresca criterion [11] 
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Figure 2.5: Graph of uniform traction bounds for a square plate with a 

central circular hole using Mises criterion [12] 




Figure 2.6: Stress (a) and displacement (b) fields for a square plate 

with a central circular hole (uniform displacement) 
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3. ELASTIC -PLASTIC BUCKLING AND POST-BUCKLING ANALYSIS 



3 . 1 Introduction 

Elastic -plastic buckling is a typical failure pattern of thin flat 
plates under in-plane compressive loading. The ultimate strength of 
structural members defines their ability to continue providing support 
after suffering damage, and is reached after compressive loading has 
caused stresses above the yield strength. The effective width concept 
for unperforated flat plates describes the limiting effect of 
slenderness upon ultimate strength. The buckling of flat plates 
weakened by circular holes has been fairly extensively studied 
theoretically and experimentally. Kurshin and Matveev [15] made a 
comparative analysis of five theoretical and three experimental studies 
of the stability of a square plate with a central circular hole. The 
best agreement between theoretical analysis and experimental results was 
achieved by Ritchie and Rhodes [16]. Ritchie and Rhodes examined square 
and rectangular plates, with simply supported unloaded edges, under 
uniform displacement compressive loading for both buckling and 
post-buckling conditions. The effective width concept discussed below 
supports their conclusion that post-buckling strength is not 
significantly affected by small or medium sized holes. 

3 . 2 Effective Width Concept 

For a plate subjected to uniaxial compressive loading, in general, 
only the parts of the plate close to the unloaded edges experience 
plastic yield. The stress distribution shown in Figure 3.1 is due to 
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the "bowing” effect in the central zone of the plate. In 1932, von 

Kdrmdn postulated that the maximum load a plate with simply supported 
edges can sustain is reached when the edge stress equals the yield 
strength [17]. The effective width of the plate b , also shown in 

e 

Figure 3.1, can be described as; 



be 

b 



_1 

Ap 



for Ap > 1 



(3.1) 



where Ap is the plate reference slenderness: 



, „ b 1.05 

Ap = u “ZI 

h /IT 



r o ^ 



1/2 



(3.2) 



and where h is the plate thickness, and K is the plate buckling 
coefficient : 



K 



f a mb 
[mb a 



(3.3) 



where a and b are the plate dimensions, and m is a wavelength parameter 
chosen to yield the smallest eignenvalue. A modern warship's strucural 
plating, disregarding the effect of stiffeners, has a typical 
length- to- thickness ratio (b/h) on the order of 200 to 300. With a 
reference slenderness ranging from 4 to 8 , such plating is relatively 
very thin. By combining Eqs . 3.1 and 3.2, the effective width can be 
described as: 



v~r r 

1.05-1 



(3.4) 



With simplified dimensions of 8 feet by 8 feet, and a thickness of half 
an inch, warship structural plating has an effective width of roughly 
1.6 feet. For plates with high slenderness values, the effective width 
is small and damage in the central zone has little impact on residual 
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strength . 



3 . 3 Buckling and Post- Buckling Analysis 

It is possible to analyze the initial instability of a flat plate 
with a central circular hole, according to Ritchie and Rhodes [16], as 
simply that of a flat plate. The strain energy U stored in a plate can 
be written as the sum of the bending and membrane strain energies: 



U - U + U 

b n 



(3.5) 



where U is the bending strain energy and U is the membrane strain 

b m 



energy . 



These strain energies can be expressed in terms of the 



out-of-plane deflection w and the membrane stress resultants N , N and 

x y 

N as : 

xy 



Ub -§ rrlfc + . 2(1-1/) [— — - U 

2 ^ [Ux 2 By 2 ' Ux 2 ay 2 Uv 

Un> - ^ Tf N £ + N £ + N £ 

zJJ I xx xx yyyy xy xyl 



w 



dxdy 



L dxdy J 
dxdy 



(3.6a) 



(3.6b) 



where D is the bending rigidity of the plate, w is out-of -plane 
displacement, and v is Poisson's ratio. Ritchie and Rhodes obtained 
the stress resultants N , N and N in each element with respect to a 

xx yy xy 

unit applied displacement of the loaded edges. The magnitude of all 
stress resultants were then related to the loaded-edge displacement by a 
magnification factor, and the strains related to the stresses by Hooke's 
law. The deflections were represented by a single sine wave in the 
loaded direction, which is a reasonable approximation for plates with 
small holes whose behavior is similar to that of unperforated plates, 
but may be less accurate for plates with large holes. 
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Ritchie and Rhodes rewrote the strain energy in matrix form in 
order to solve for the buckling modes and loads. The non- trivial 
solutions of the simultaneous equations minimizing the strain energy 
were the eigenvalues of the edge displacement coefficient <j> and the 
corresponding eigenvectors. The solutions also satisfy the criterion 
for neutral equilibrium. Thus the smallest eigenvalue and its 
corresponding eigenvector give the buckling load or displacement and the 
buckling mode. Ritchie and Rhodes used a finite element analysis, with 
a mesh of 97 triangular constant- strain elements over one quarter of the 
plate, to obtain the stress distribution for both uniform edge 
displacement and uniform edge stress cases. The uniform displacement 
and uniform traction loading curves thus obtained are presented in 
Figure 3.2. 

The post-buckling behavior was examined by a combined semi -energy 
and finite-element method. It was assumed that post-buckling 
out-of-plane deflections remain of the same form, but with varying 
deflection magnitude, as those obtained from the buckling solution. 
This assumption was considered reasonable for plates without holes if 
the applied load was less than about twice the buckling load. The 
presence of local buckles, however, causes changes in the stress system 
throughout the plate. These changes can be evaluated from the solution 
of von Kdrmdn's comparability equation, which relates in-plane and 
out-of-plane deformations in a plate by means of a stress function F, as 
follows : 
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V F 



rr_aV ) 2 a 2 w a 2 w ] 
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(3.7) 
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A second system of stresses was derived by solving von Kdrmdn's 
comparability equation, and applying the boundary conditions at the 
plate supported edges. A third stress system was obtained, using finite 
elements, by hypothetically applying a loading system around the hole 
edge to counteract the stresses produced by the second system. The 
complete membrane stress system in the plate was obtained from the 
superposition of these three systems: 

n =<£n +C 2 fn +n *) (3.8a) 

xx xxl r ^ xx2 xx3 J 

n - <f> n +C 2 fn + n ) (3.8b) 

yy yyl r ^ yy2 yy3 J 

n -= <£ n + C 2 f n + n 1 . (3.8c) 

xy xyl r ^ xy2 xy3 J 

where the subscripts 1 , 2 , and 3 refer to the respective stress system. 

The principle of minimum potential energy was applied by Ritchie 
and Rhodes to their post-buckling analysis in order to determine the 
buckled form. A change in natural wavelength of the rectangular plate 
buckled form, from two to three half -wavelengths , occurs if the hole 
diameter is greater than 0.35 times the plate width. Brush and Almroth 
[18] developed a plot of buckling coefficient K versus aspect ratio 
(a/b) , as shown in Figure 3.3, which presents eigenvalue curves for 
critical values of axial stress for simply supported unperforated plates 
subjected to in-plane compressive loading. A simply supported square 
plate would be expected to have a buckled form of one half -wavelength . 

3 . 4 Ultimate Strength of A Holed Plate 

A structure can continue to support loading after buckling. As 
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with von Kdrmdn's analysis, the failure of a plate under a compressive 
edge load is assumed to occur when the maximum stress in the loaded 
direction reaches the yield stregth. For a simple analysis, Ritchie and 
Rhodes defined the ultimate load as: 



P - a b h (3.9) 

ult Ye 

where b is the effective width corresponding to the yield stress a^, 
which is similar to Eq . (2.5). The effect of a hole in a buckled plate 

on the residual strength can be approximated by redis tr ibuting the 
compressive stresses released by the introduction of the hole over the 
remaining region. This produces an increase in maximum stress for a 

given load, and a corresponding reduction in effective width. The 
post -buckling stress distribution can be represented by a simple sine 
wave as shown in Figure 3.4 (which is similar to that shown in Figure 
3.2). The increase in maximum stress due to the hole is calculated by 
redistributing the compressive stresses "removed" by the hole: 

P «= h f a dy (3.10) 

h J x 
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where c is the region carrying tensile stresses 

plate, c is the diameter of the hole, and P is 

h 

removed by the hole. The change in stress is: 
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For a given external load P: 



P — o b h = (a + A a ) b h (3.12) 

xmaxe x x ec 

where (a + A a ) is the maximum stress after the introduction of the 

X X 

hole and b is the corrected effective width. The corrected effective 
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width can then be written as: 



b 

e c 



b 




(3.13) 



Figure 3.5 was produced using the relationship between applied load and 
effective width separately derived by Rhodes and Harvey [19]. The 
ultimate load for plates with various size holes was then calculated 
and plotted in Figure 3.6, which indicates that ultimate strength is 
essentially unaffected by the presence of small holes (diameter less 
than half the width) and is roughly equal to three times the critical 
buckling load P . 

No attempt was made to redistribute tension stresses released by 
the introduction of the hole into the analysis. This is a conservative 
approximation because the effect of redistributing high tensile stresses 
could indicate that the effective width of a plate with holing is 
greater than that of an unperforated plate. While it is possible that 
this could indeed be the case in practice, Ritchie and Rhodes suggest 
that it would be dangerous to advocate such a hypothesis for use in 
design without experimental proof. 

Roberts and Azizian [9] have examined this question further and 
state that for uniaxial and biaxial compression the buckling load may 
increase with increasing size of the central circular hole due to the 
redistribution of membrane stresses towards the edges of the plate. The 
increase in buckling load is dependent upon the boundary conditions and 
appears significant only for plates with clamped edges. Roberts and 
Azizian conclude that "even though the buckling load may increase with 
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increasing size of the hole, it is not expected that there will be a 
corresponding increase in the ultimate or collapse load." 



3 . 5 Comparison With Experimental Results 

The experimental investigation of Ritchie and Rhodes [16] consisted 
of a series of tests on square and rectangular hot-dip galvanized mild 
steel plates to collect information on buckling loads, out-of -plane 
deflections, collapse loads, and strain distributions along the axes of 
the plates. Plate dimensions were nominally 254mm wide by 254mm or 
508mm long (square and rectangular respectively) and 1.57mm thick. The 
unloaded vertical edges of the plate were supported in the frame by 
knife edges and loaded edges were supported in rollers resting on needle 
bearings contained in the blocks, through which the load was 
transmitted. The specimens were thus simply supported on all edges and 
subjected to uniform displacement loading. The buckling loads were 
estimated by obtaining the intersection of the pre- and pos t-buckling 
slopes of the load membrane strain plots. All plate specimens were 
loaded under uniform displacement conditions. 

Ritchie and Rhodes normalized their analytic solutions and 
experimental results with respect to the critical buckling load. From 
Timoshenko and Gere [20], the critical stress of a square plate with 
similar boundary conditions can be described as: 



2.528 h E 

T = 

crit 2 N 2 

(1 - " ) a 



( 3 . 14 ) 



The critical buckling load can thus be normalized with respect to the 
tensile yield load as: 
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(3.15) 
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For the mild steel square plates investigated by Ritchie and Rhodes, 
assuming a yield strength of 30 ksi, the normalized critical buckling 
load is roughly 0.106. 

The analysis developed by Ritchie and Rhodes describes a normalized 
ultimate load of roughly 0.30, with respect to the tensile yield load, 
for their square plates. Both the critical stress and the effective 
width are strongly affected by relative plate thickness. The effective 
width concept describes a normalized ultimate load for an unperforated 
square plate, of the same dimensions and material properties as tested 
by Ritchie and Rhodes, of roughly 0.34 with respect to the tensile yield 
load. 



The theoretical relationship between buckling load and hole size in 
simply supported plates loaded by both uniform traction and uniform 
displacement is supported by the experimental results. The experimental 
results of other researchers for uniform traction buckling of a square 
plate are plotted in Figure 3.7 along with the analytic solution 
developed by Ritchie and Rhodes. The experimental results for uniform 
displacement buckling of a square plate are plotted in Figure 3.8, along 
with the analytic solutions. The change in buckling mode in rectangular 
plates from two to three half-waves, when the hole diameter exceeds 0.35 
times the plate width, was observed in the experimental investigation. 
The scatter in experimental results, according to Ritchie and Rhodes, 
was due to the large initial imperfections in the plates. The ultimate 
uniform displacement loads for both square and rectangular plates are 



38 



presented in Figure 3.6 and indicate good correlation between the 
analytic solution and experimental results. The ultimate strength of a 
thin plate is essentially unaffected by the presence of small holes, and 
is roughly equal to three times the critical buckling load. 

With regards to application, however, the experimental results 
indicate that the elastic -plastic ultimate strength analysis developed 
by Ritchie and Rhodes is slightly conservative. The observed normalized 
ultimate load for a thin unperforated square plate ranged from 0.34 to 
0.35, which is most accurately described by applying the effective width 
concept. It is suggested that the Ritchie and Rhodes analysis, which 
accurately predicts the loss of strength associated with increased hole 
size, be "shifted" to coincide with the effective width solution for an 
unperforated plate. 
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Figure 3.1: Real stress distribution and effective plate width due to 

in-plane compressive loading [18] 




Figure 3.2: Graph of theoretical uniform traction (constant stress) 
and uniform displacement (strain buckling) buckling loads versus cutout 
ratio for square plates [16] 
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Figure 3.3: Buckled form and critical values of axial stress for simply 

supported plates subjected to in-plane compressive loading [18] 





Figure 3.4: Assumed stress distribution across holed plate [16] 
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Figure 3.5: Graph of effective width versus maximum stress for various 

cutout ratios [16] 




Figure 3.6: Graph of theoretical analysis and experimental results for 

ultimate load versus cutout ratio for square and rectangular plates [16] 
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Figure 3.7: Graph of theoretical analysis and experimental results for 
uniform traction (constant stress) buckling load versus cutout ratio for 
a square plate [16] 
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Figure 3.8: Graph of theoretical analysis and experimental results for 
uniform displacement (constant strain) buckling load versus cutout ratio 
for a square plate [16] 
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4. DISHING OF PLATES UNDER IMPULSIVE LOADING 



4 . 1 Introduction 

Shipboard structural members are deformed by the over- 
pressurization associated with an internal explosion. The effects of 
holing produced by the initial shock wave upon residual strength of a 
flat plate have been investigated in Chapters 2 and 3. The quasi -static 
overpressurization which follows the blast produces lateral deflections 
in the bulkheads and deckplates of compartments . 

The deformation of thin flat plates subjected to large impulsive 
loading has been studied extensively since World War Two. Thorough 
reviews of the historical development of empirical and analytical 
studies have been presented by several authors [21,22,23]. A number of 
researchers have developed good analytic techniques. In structural 
elements such as beams and plates, bending resistance increases along 
with membrane forces while the deflection is much less than the 
thickness of the element. As the deflection increases, however, 
stresses due to stretching of the middle plane dominate over bending 
stresses. At this stage the behavior of a beam is comparable to that 
of a plastic string. Perrone and Bhadra [24] developed this concept 
with a simple lumped mass at the midpoint. Recently, Nurick [23] 
developed a highly accurate membrane mode analysis based upon the lumped 
mass string concept. Excellent empirical relationships between 
deformation and impulse are developed from extensive testing conducted 
by Nurick and others. 
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4 . 2 Dynamic Effects of an Internal Airblas t 

Roughly 70 percent of the energy released by a warhead detonation 
produces blast effects. The critical loading associated with an 
explosion is the high pressure discontinuity at the shock front. The 
shock wave produces deflection, necking, rupture and holing of plating 
close to the point of detonation due to the high velocity gradients 
across the plate surface. Rupture begins at the point of a stiffened 
plate structure nearest the explosion and can extend to the plates' 
boundaries. The shock loading is predominately normal to all structural 
members . 

After passage and decay of the shock wave from an internal 
airblast, confinement of the explosion gases leads to a pressure loading 
which has a long duration compared to the response time of typical 
stiffened structures. There is no simple relation between explosive 
mass and pressure impulse [23]. As a result of damage to bulkheads, 
decks, and overheads surrounding the burst space, the gas overpressure 
is vented to the atmosphere and adjacent spaces and is relieved as the 
confinement volume increases. This "quasi - static” pressure loading is 
governed by the confinement volume. Figure 4.1 shows a typical trace of 
the pressure versus time history recorded at a point on the wall of a 
partially closed chamber following an internal explosion. Pressure 
peaks are caused by repeated reflections of the blast wave off the 
chamber wall. The impulse imparted by the pressure peaks is small 
because of their short duration, and produces negligible damage effects 
on internal structure and equipment [6,25]. 
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4 . 3 Membrane Mode Finite Element Analysis 



Membrane mode approximations, in which bending is ignored, have 
been applied to impulsively loaded thin plates by several authors who 
all assumed a fixed mode shape. Nurick [23], however, has developed a 
model which considers lateral displacements in addition to transverse 
deflections, and which computes the mode shape at each time step. 

The string concept forms the basis of the Nurick model. The 
string is sub-divided into several sections with associated lumped 
masses, and the acceleration is assumed to be a linear function of 
the velocity. As developed by Perrone and Bhadra [24] , the 
acceleration is the second differential of displacement resulting in a 
first order nonlinear differential equation which is solved by numerical 
integration. Transverse and lateral deflections are incorporated 
into the equation such that the trajectory of every point on the 
surface remains normal to the plate surface. 

The material is assumed to behave in a rigid visco-plastic 
manner. This can be described by the empirical relationship developed 
by Bodner and Symonds [26]: 




where e is the strain rate, D and n are material constants. Manjoine 
[27] investigated the influence of temperature and rate of strain on the 
yielding of mild steel. Bodner and Symonds applied Manjoine 's work to 
their empirical relationship in order to define the material constants. 
For mild steel the material constants D and n are 40 s ^ and 5 
respectively . 

The deformed geometry of the string can be defined by a modeshape 
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<f ) , which varies with space and time. The transformed displacement is 
written as: 

w (r , t) = 4> (r , t) w (t) (4.2) 

i l c 

where (r,t) are the coordinates of space and time and w (t) is the 

c 

central displacement. It is further assumed that the transverse 
velocity and acceleration fields at time t can be differentiated 
without taking into account the variation of the mode shape <f> with time, 
which is known as an instantaneous or stationary mode assumption. 

The string concept is applied to plates by discretizing the plate 
into contours along which <f> is constant. The mode shape 4> is assumed to 
vary linearly with distance between contour levels. Symmetry dictates 
square contours for the square plates, as shown in Figure 4.2 . This 

symmetric contour is justified by the deformation shape, shown in Figure 
4.3, produced by explosive impulse loading. 

It is possible to write (t + 1) dynamic equations for the lumped 
masses. Nurick developed an algorithm to solve the dynamic equations. 
The initial central velocity is computed from the initial total impulse 
I, which is assumed to impart a uniform velocity to the plate. The 
initial mode velocity is then obtained from a generalised momentum 
balance. The stiffness of the plate is assumed to be zero for the first 
time step, and the iterative solution commences with the second time 
step. The equations are highly non-linear, and are solved using a 
Newton-Raphson technique. The iterative procedure continues until a 
prespecified tolerance is met in order to continue to the next time 
step, and the forward integration is terminated when w^ equals zero. 
Convergeance of the algorithm suggests that the final shape is reached 
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within 100/is. 



4.4 Experimental Results and Empirical Analysis 

The experimental investigation of Nurick [23] consisted of a series 
of tests on circular, square and rectangular plates to collect a range 
of deflection information on plates subjected to high impulsive loading. 
Plate specimens were cut from 1.6mm thick cold rolled steel and were not 
annealed. Plate dimensions were nominally 100mm diameter circular, 89 
by 89mm square, and 113 by 70mm rectangular. Metabel was used as the 
sheet explosive for the experiments, and was placed on a 16mm thick 
polystyrene pad to attenuate the shock transmitted to the plate, provide 
a uniform impulse, and prevent spallation of the specimen. The specimen 
plate was clamped between two holding plates attached to a ballistic 
pendulum used to measure the impulse. A light interference method 
was used to measure the deflection- time history. The final 
mid-point deflection was measured by means of a vernier and also by a 
reflex metrograph. 

The range of experimentation performed by Nurick was impressive. 
Experimentation by previous researchers did not extend beyond 
deflections of six to seven plate thicknesses for circular plates and 
eight to nine thicknesses for quadrangular plates. Nurick performed a 
total of 140 impulse loading experiments with deflections ranging up to 
twelve plate thicknesses. 

The analytic prediction method developed by Nurick compares 
favorably with impulse loading experimental results for circular, 
square, and rectangular plates. The range of analytic predictions and 
empirical results is shown in Figure 4.4 for circular plates and Figure 
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4.5 for square plates. Although most of the prediction methods can do 
an adequate job predicting deflections with small deflection- to- 
thickness ratios, they display tendencies to highly over or under 
estimate deflections with large ratios. These over and under estimates 
are based upon predictions which assume a fixed mode shape. 

For circular plates, the best agreement between predictions and 
experimental results was achieved by Nurick, Lippman [28], Perrone and 
Bhadra [29], and Symonds and Wierzbicki [22]. The closed- form solution 
for circular membranes, developed by Symonds and Wierzbicki, is a simple 
but accurate approximation method which disregards bending. For 
impulses which blanket the plate surface, their solution can be 
simplified as: 



w 



c f 



.212 I 

■ h 0 O 
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(4.3) 



where w is the final deflection at the center of the plate, r is the 
cf 

plate radius, h is the plate thickness, and p is the material density. 
For quadrangular plates, the best agreement between predictions and 
experimental results was achieved by Nurick, and Baker [30]. Using an 
energy method, Baker assumed that the energy imparted to the structure 
could be equated to plastic work which he then equated to the initial 
kinetic energy for impulsively loaded structures. Baker assumed that 
the force or pressure pulse was applied in a time short enough for the 
structure response to be only a function of the integral of the pulse, 
rather than the details of its time history. 

It is interesting to compare and, if possible, utilize the 
historical data available from different impulse loading 
experiments. A dimensionless damage number was defined by Johnson [31] 
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as : 



a 



2 

£v_ 

a 



(A. 4) 



o 

where v the impact velocity and o is the damage stress (mean flow 

o 

stress) . Table 4.1 shows the regime of damage as a function or damage 
number, and shows the results of some circular plate experiments. The 
Johnson damage number was modified by Nurick in order to 

non-dimensionalize dissimilar geometries. A modified damage number for 
circular plates can be written as: 
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and a modified damage number for quadrangular plates can be written as: 
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where a and b are breadth and length respectively of a quadrangular 

plate, and a is given the value of the static yield strength a . 

o y 



Table 4.1: Typical Johnson's Number Damage Regime 

(mild steel at room temperature impacted by a projectile) [31] 
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For plates of equal area and thickness subjected to equal total 
impulses , the midpoint deflection of the circular plate exceeds 
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that 



of the 



that of the square plate, which in turn exceeds 

rectangular plate. A dimensionless plot of deflection- to - thickness 

ratio versus modified damage number 6 for quadrangular plates is 

q 

presented in Figure 4.6. A least squares correlation analysis of 
the experimental results indicates that Nurick's experimental results 
are consistent with previous data, and that a relationship exists 
between the impulse a given plate is exposed to and the deformation it 
experiences . 

Based upon the extensive range of consistent experimental 
results, it is possible to develop an empirical solution to the 
deflection problem. A least squares analysis of the circular plate 
experimental data yields: 

T 1 = 0.425 0 c + 0.277 (4.6a) 



where the number of data points was 109 with a correlation coefficient 

of 0.974. A least squares analysis of the quadrangular plate 

experimental data yields: 

w 

- 0.236 6 + 0.001 (4.6b) 

h J q q 

where the number of data points was 156 with a correlation coefficient 
of 0.984. The empirical relationships can be simplified to: 

- 0.425 0 (4.7a) 

c c 

= 0.236 B (4.7b) 

q q 

as a useful "engineering" approximation. 

The excellent correlation of the empirical formulas indicate that 

they can be readily applied to engineering analyses. The linear nature 

of the closed-form deformation solution developed by Symonds and 




51 



Wierzbicki, and others, supports the use of a linear empirical 
relationship based upon a modified Johnson's damage number. Extensive 
deformations can produce tearing along the edges of the plate or holing 
in the center of the plate. The maximum non- tearing 

deflection- to- thickness ratios during Nurick's testing were 11.60 and 
11.97 for circular and quadrangular plates respectively. This tearing 
along the clamped boundaries indicates that both analytical and 
empirical prediction methods are limited to modified damage values of 6 

c 

< 25 for circular plates and 6 <50 for quadrangular plates. 
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Figure 4.1: Graph of pressure versus time response recorded at a point 
on the wall of a partially closed chamber following an internal 
explosion [25] 




Figure 4.2 Analytic model showing discretization of a square plate [23] 
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Figure 4.3 Final mode profile and contour plot of a typical explosive 
test for a square plate [23] 




Figure 4.4: Graph of analytic solutions and experimental results for 

deflection- to- thickness ratio versus impulse for circular plates [23] 
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Figure 4.5: Graph of analytic solutions and experimental results for 

deflection- to - thickness ratio versus impulse for square plates [23] 




Figure 4.6: Graph of deflection- to- thickness ratio versus dimens ionless 

damage number 6 for quadrangular plates [23] 
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5. STRENGTH OF A DISHED PLATE 



5 . 1 Introduction 

Despite the large lateral deformations suffered by shipboard 
bulkheads and deckplates as a result of impulsive loading, the plating 
can continue to support normal service loads. The dishing of thin flat 
plates into a shell was related to such parameters as material 
properties and explosive impulse in Chapter A. By examining large scale 
elastic -plastic deformation mechanisms it is possible to describe the 
buckling behavior of a thin shell. 

An out-of- flatness plate with a small initial deflection does not 
react to in-plane loading like an ideal flat plate. The effect of 
initial out-of-plane deformations is generally to reduce the peak load 
of the plate and bias the mode of failure towards the initial shell 
shape. The effect of initial imperfections upon the ultimate strength 
of plates offers some insight into the effect of large scale dishing. 
The residual strength of a thin shell with large initial deformation (on 
the order of 5 to 15 plate thicknesses) is first analyzed as a simple 
extension of the concept of initial imperfections. An extensional 
collapse mode analysis is then developed as a simple closed-form 
solution which accounts for actual deformation mechanisms. 

5 . 2 Physical Basis of Deformation Mechanisms 

The yield load of a plate and the elastic -plastic buckling load 
provide two limits to plate strength. The onset of plasticity is 
affected by both the in-plane and plate bending stress components. 
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Plasticity becomes more significant as the initial out-of-plane 
imperfections become larger causing surface and membrane yield to occur 
earlier. The interaction of membrane and bending stresses, with initial 
imperfections, produces an early loss of stiffness and reduction in peak 
load below the ideally planar elastic buckling load [17]. 

Although the ultimate strength of a plate with very slight initial 
imperfections remains reasonably close to that of an ideal flat plate, 
the effects of larger imperfections (on the order of the plate 
thickness) can be significant [17]. Von Karman developed the effective 
width concept (presented in Chapter 3) to describe the stress 
distribution produced in a simply supported plate subjected to uniaxial 
compressive loading. The influence of geometrical imperfections has 
been incorporated into a von Karmdn analysis by several authors, as 
shown in Table 5.1. Initial deflections of the plate lead to a marked 
reduction of the peak load, especially for thin plates with reference 
slenderness Ap less than 1.5. From Eq . 3.2, the reference slenderness 



is : 



„ b 1 - 05 
p = 




(5.1) 



Rectangular plates restrict the failure mode and generally have greater 
stiffness than square plates, but suffer a steeper post-buckling 
behavior [17]. 

In order to understand the true collapse behavior of a shell 
subjected to uniaxial compressive loading, the effect of out-of-plane 
deformations on the in-plane stress state of the plate must be 
considered. In the large deflection regime, in general, in-plane and 
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out-of -plane displacements and hence stresses and strains interact. As 
a result of dishing, the length of the plate center line becomes greater 
than the boundary length and membrane tensions develop in both plate 
directions. Membrane action becomes more and more dominant as 
out-of-plane deflections become larger [17]. 

Table 5.1: Empirical relationships to describe the effect of 
small initial imperfections upon effective plate width [17] 
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5 . 3 Large Initial Imperfection Analysis 

Extensive research has been conducted in the last twenty- five years 
on the ultimate load behavior of box girder structures [32,33,34]. This 
research has followed a series of structural failures involving merchant 
ships (large ore carriers in rough seas near Japan) [33] and large box 
girder bridges in Europe [34]. Frieze, Dowling and Hobbs [32] have 
developed a large deflection elastic-plastic analysis to examine the 
behavior of plates in compression. 

The analytic method developed by Frieze et al . uses a modified von 
Karman analysis to account for material non-linearity and the effects of 
residual stresses and initial out-of-plane deflections. The governing 
equations were solved by computer processing using an iterative finite 
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difference procedure . 



Frieze et al . performed a parametric study 



examining the collapse behavior of a uniaxially loaded plate. This 
analysis indicates that the influence of aspect ratio (a/b) is dependent 
upon both the magnitude of the initial deformation and the slenderness 
ratio. Although the in-plane restraint along unloaded edges 
(constrained or unrestrained) has little effect on the behavior of 
stocky plates, the analysis indicates that slender plates with 
constrained edges retain their strength better than with unrestrained 
edges. Figure 5.1 shows that this is especially true for plates with 
large initial deformations. Slender plates, in particular, suffer from 
the weakening effects of initial deformations which generally dominate 
the effects of residual stresses. 

The residual strength of a thin shell can be extrapolated from 
theoretical and empirical analyses which have been performed for modest 
initial deformations. It is not possible to directly apply the analyses 
because shipboard structural plating is very thin. The analyses can be 
extrapolated by applying the von Kdrmdn effective width concept to 
describe the loss of strength for unperforated flat plates with 
increasing slenderness. The theoretical analysis of Frieze et al. is 
summarized by Figure 5.2, which presents six different curves of 
ultimate load versus slenderness ft for different levels of initial 
dishing. The slenderness value f) is defined as: 



^ h -I E 



a Y 



(5.2) 



which is roughly twice the plate reference slenderness Ap for a square 
plate. The slenderness ratio for shipboard structural plating, 
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disregarding the effect of stiffeners, ranges roughly from 9 to 15. 
Ueda et al. [35] also examined the effects of initial deformations and 
residual stresses on the compressive strength of square plates, and 
conducted a modest program of experimentation. Large-deflection, 
elastic -plastic analyses were performed using an incremental finite 
element technique. Ueda's theoretical and experimental research 
supports Freize's analysis, as shown in Figure 5.3, and by 
extrapolation suggests that the normalized ultimate load would decrease 
for a thin shell with large initial deflections by roughly 10 percent 
from that of a flat plate with similar slenderness. The plate critical 
buckling curve is included in both Figures 5.2 and 5.3 for comparison 
purposes . Even with large initial imperfections, the ultimate strength 
of slender plates is greater than the critical buckling load. 

5 . 4 Extensional Collapse Mode Analysis 

A structure can absorb crushing work through local plastic folding 
with bending and extensional deformation. Minorsky [36] examined the 
ship collision problem, with an analysis of empirical data from several 
dozen reported casualties, and found that plastic work dissipation is 
linearly related to the volume of damaged members. Although useful for 
studying major collisions, Minorsky' s method is not valid for structural 
failure of a local nature. The theory of extensional collapse of plate 
elements was developed by Hayduk and Wierzbicki [37] and others in the 
early 1980 's to examine the crush behavior of individual members during 
crushing. Extensional collapse mode analysis is described here in order 
to determine the average load-carrying capability in structural plating 
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crushed beyond the ultimate load. 

Ships are generally constructed of plate-like structures, as shown 
in Figure 5.4. During deformation, the edges of the structure can be 
considered hinge lines, as can the horizontal limits of the buckling 
region. Crushing work is dissipated through plastic bending of adjacent 
plates about the horizontal hinge lines connecting them and through 
shape distortion of each plate. The energy-balance postulated by 
Wierzbicki and Abramowicz [38] indicates that zones of bending and 
extensional deformation, which are restricted to a small fraction of the 
total area of the structure, dissipate essentially all of the work done 
by the applied load during the crushing deformation mode. 

In the propagating hinge line deformation mode (Mode I) , extension 
is limited to the neighborhood of the hinge line. Displacement fields 
must be described by a continuous and smooth function across a 
propagating hinge line. The quasi- inextensional collapse mode of a re- 
shaped angle element is shown in Figure 5.5. Wierzbicki and Abramowicz 
described the individual dissipation energies associated with the 
deformation of a toroidal shell, horizontal hinge line, and inclined 
hinge line through complete folding of the section. In the basic 
folding mechanism of a L- shaped angle element there are eight horizontal 
hinge lines, two inclined hinge lines, and one toroidal shell section. 
The resulting theoretical crushing strength of a L-shaped element is 
thus given by: 




where P is the mean crushing load and M is the fully plastic yield 
m o 
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moment defined as: 




2 



(5.4) 



where the flow stress o can be reasonably approximated as: 



o -=0.9 o 



UTS 



(5.5) 



O 



Short- angle columns commonly fold up and fail in an asymmetric mode 
unless very large imperfections are introduced in the shape of the 
symmetric modes [37]. 

In the stationary hinge line deformation mode (Mode II) the edges 
of the structure remain stationary while the side plates undergo 
rigid-plastic distortion. Plate distortions consist of either extensive 
in-plane compression or in-plane extension. Hayduk and Wierzbicki 
assumed that there is a linear variation of the velocity field in the 
horizontal y-direction, as shown in Figure 5.6, and developed the 
individual dissipation energies associated with the deformation of a 
plate, stationary horizontal hinge and stationary inclined hinge through 
complete folding of the section. 

Hayduk and Wierzbicki applied the extensional collapse mode 
analysis to a x- shaped crucifom structural element. A cruciform is 
constructed of four flanges which are rigidly joined along the initially 
vertical edge. A cruciform element reasonably simulates the boundary 
conditions on a plate in a typical shipboard compartment. Within the 
folding zone, the joint line becomes an inclined hinge fixed in the 
material. In Mode II crushing of a cruciform there are eight 
trapezoidal distortion sections, sixteen stationary horizontal hinge 
lines, and four stationary inclined hinge lines. An upper bound to the 
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crushing strength of a cruciform element can thus be described as: 



M 



P 



m 



o 



- 20.05 




(5.6) 



This least upper bound is based upon a minimum energy analysis which 
indicates that folding will occur with a zone of half -wavelength H: 



The upper bound analysis also indicates that the energy dissipated at 
the inclined hinge lines is negligible for a sufficiently large 
length- to- thickness ratio . 

A lower bound on the crushing strength of a cruciform can be 
similarly derived by noting that a continuous joint along the vertical 
edge of a cruciform increases the column strength as compared to a 
column composed of two unjoined L-shaped elements. From the Mode I 
analysis of a single L-shaped element, a lower bound to the crushing 
strength of a cruciform can be described as: 



Careful inspections of several crushed cruciforms, such as the one shown 
in Figure 5.7, reveal that the actual collapse mechanism is a 
combination of the two basic deformation modes. Hayduk and Wierzbicki 
developed a modified Mode II model, and described the mixed-mode 
crushing strength of a cruciform as: 



The Hayduk and Wierzbicki extensional collapse mode solutions for 




(5.7) 
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(5.8) 




(5.9) 
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cruciform failure are based upon complete folding of the structure. 

This extensional collapse mode analysis of a cruciform can be 
applied to the residual strength of a dished plate. The typical 
load-deflection characteristics of a crushed cruciform, shown in Figure 
5.8, indicate that the ultimate strength of a cruciform with initially 
flat side plates is rapidly reached during uniaxial loading. The 
post-buckling and folding behavior after the ultimate strength is 
reached "reinforces" the cruciform and allows it to continue to support 
loading during repeated folding. Since the collapse mode analyis 
indicates that folding will occur along a zone of width 2H defined by 
Eq. 5.5, dishing may not significantly reduce a plate's ultimate 
strength because the dished shape does not encourage the folding. 
Although dishing may reduce a plate's ultimate strength, it is 
postulated that the mean crushing load P which the structure can 

m 

support is little affected by the dishing. 

The extensional collapse mode analysis equations can be rewritten 
to provide normalized loads for shell crushing. The mean crushing load 
is normalized with respect to the tensile yield load. By assuming that 
the strength of a single dished plate can be considered equal to 
one-half that of a cruciform, the upper bound can be expressed as: 

P a 2.51 h (b/2h) 1/2 



P a 

Y Y 



(5.10a) 



the lower bound can be expressed as: 



P cr 3.01 h (b/2h) 



1 / 3 



P a 

Y Y 



(5.10b) 



and the mixed-mode solution can be expressed as: 
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P o 2.87 h [b/2h] 1/3 + 0.36 h (b/2h) 



2 / 3 



P O 

Y Y 



(5.10c) 



A modern warship's structural plating, disregarding the effect of 
stiffeners, has a typical length- to- thickness ratio (b/h) on the order 
of 200 to 300, From the upper and lower bound analysis for a warship's 
structural plating, the normalized mean crushing load ranges from 0.05 
to 0.13. 



The extensional collapse mode analysis can be applied to a holed 
shell. Since the width of the local folding zone is much smaller than 
the height of the typical cruciform, as shown in Figure 5.9, the "wing 
areas" can be considered ineffective in terms of both extensional and 
bending energy. By assuming that the fold is made through the hole 
region, the effective plate half-width can be reduced by the radius of 
holing. For Eqs. 5.10a - 5.10c only the half -width- to - thickness ratio 
(b/2h) in the numerator is reduced by holing. The extensional collapse 
mode analysis can then be applied using the reduced half -widths. 



5 . 5 Comparison With Experimental Results 

The experimental investigation of Hayduk and Wierzbicki [37] 
consisted of a series of tests on thin copper cruciform specimens. 
Specimens were constructed of two L-angle elements of 0.014 inch thick 
copper sheeting butted along the bend line and soldered with common 
60-40 lead-tin solder. Specimen length varied from 2.0 to 4.0 in. while 
plate half-widths varied from 1.0 to 2.0 in., for a length- to - thickness 
ratio which varied from 143 to 286. The unloaded vertical edges of the 
cruciform were clamped at the top and bottom to prevent specimen 



65 



twisting during loading- During quasi-static loading the specimens were 
routinely crushed from 40 to 60 percent of their original length. 

The extensional crushing mode analysis is supported by the 
experimental results. Figure 5.10 presents the experimental data and 
the upper, lower and mixed-mode solutions for crushing of copper 
cruciforms and other elements. The experimental data for cruciforms are 
shown by open circles. The scatter of experimental data indicates that 
the cruciforms failed according to a varying combination of collapse 
mechanisms . The collapse mechanisms are sensitive to imperfections . 
Indications of fracture at the joint in some test pieces explains some 
of the scatter below the mixed-mode solution. Despite the scatter of 
data, all experimental results for cruciforms lie within the upper and 
lower bound solutions to normalized mean crushing load. Amdahl [39] 
performed a similar series of tests on larger steel cruciforms which 
also support the validity and usefulness of extensional cushing mode 
analysis . 
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5.1: Graphs of theoretical ultimate load versus strain for 

plates with different slenderness values [32] 




5.2: Graph of theoretical ultimate load versus slenderness for 

constrained plates with different initial deflections [32] 
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Figure 5.3: Graph of theoretical analysis and experimental results for 
ultimate load versus slenderness for square plates including the effects 
of initial deflections and welding residual stresses [33] 
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Figure 5.4: Typical structural elements 
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Figure 5.5: Basic folding mode of a L- shaped angle element with a 

continuous and smooth displacement field [37] 
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Figure 5.6: Plastic distortion of a rectangle into a trapezoid [37] 
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Figure 5.7: Photograph of partially crushed cruciform 




6 . in. 



Figure 5.8: Graph of typical load versus deflection characteris t ics of 

a crushed cruciform [37] 
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Figure 5.9: 



Cruciform with holed half-plates 




Figure 5.10: Graph of theoretical analysis and experimental results for 
mean crushing strength versus half -plate length- to- thickness ratio for 
cruciforms [37] 
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6. MIXED MODE CRACK PROPAGATION 



6 . 1 Introduction 

The hole edge produced during the fracture and tearing of metal 
plates associated with an internal airblast is not smooth. Although the 
hole can reasonably be considered to be circular for some forms of 
analysis, it is important to examine the effect of a ragged and torn 
edge upon residual strength. Initial cracks of various sizes will 
emanate from around the rough edge of the hole. Crack growth under 
loading (monotonic or cyclic) may limit the residual strength of the 
structure . 

The initiation of cracks during holing and overpressurization 
associated with the airblast is a complex dynamic situation. The 
deformed plate can be considered a thin shell. Crack initiation and 
propagation from the torn edges associated with holing and fragmentation 
occurs in a fully plastic or elastic -plastic deformation field. 

Crack propagation occurs as a mixture of three limiting modes, as 
shown in Figure 6.1. Mode I crack propagation is associated with 
in-plane tension, Mode II is associated with in-plane shear, and Mode 
III is associated with out-of -plane shear. Actual crack propagation is 
rarely so distinct (e.g. slant fracture consisting of Modes I and III 
can be caused simply by Mode I loading) . As one could imagine tearing a 
page out of a telephone book, the metal plating will be opened by 
in-plane tension and torn by out-of -plane shear. Crack propagation from 
the rough edge of a hole in a thin shell will be caused primarily by 
slant fracture. 
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6 . 2 Physical Basis of Plastic Fracture 

A fully plastic situation consists of plane strain, plane stress, 
or intermediate modes which occur when the plastic zone is large 
relative to the plate thickness. The plastic zone will, for plane 
stress, become large and spread through the whole cracked section if the 
net section stress is of the order of the yield strength. In practice, 
linear elastic fracture mechanics (LEFM) does not apply if the fracture 
stress is greater than two-thirds of yield [40]. Plane stress requires 
the vanishing of one of the principal components of stress, o , and the 
vanishing of stress gradients in the z direction. Plane stress occurs 
when the out-of-plane dimension is far thinner than that required for 
plane strain, as in a thin plate or shell. 

The tearing of thin sheets of ductile metals is a process of 
combined flow and fracture. Plasticity precedes and accompanies crack 
propagation. Necking may occur when the plastic zone size becomes 
larger than the plate thickness. Stress and strain distributions in the 
region of the tip of a sharp crack are complex. Stress changes are very 
large within one plate thickness of the crack tip, and a state of plane 
strain may be approached. Plane stress behavior will be exhibited 
beyond roughly four plate thicknesses of the crack tip [10], 

Crack propagation can be regarded as a process of hole nucleation, 
coalescence, growth, and localization. Localization may occur first in 
plane stress. Tipper [41] examined sections near the fracture surface 
cut from different mild steel specimens. In ductile steels, hole 
nucleation almost certainly occurs by the drawing away of the metal from 
non-metallic inclusions. Separation of the metal around inclusions 
occurs at some distance (1 or 2 hole spacings) ahead of the crack tip, 



73 



and may occur with stresses as low as yield. The holes become linked by 
elongation and eventually produce macrofracture. 

6 . 3 Analytic Approaches to Plastic Fracture 

Plastic fracture is a complex process whose mechanisms are not 
fully understood. These mechanisms include hole nucleation, 
coalescence, growth, and localization of plastic flow into bands. 
Linear and non-linear elastic fracture mechanics are relatively well 
understood for crack initiation in plane strain, but do not apply to the 
loading of shells and subsequent crack growth. Despite considerable 
research within the last ten years, according to Broek [40], no 
generally accepted approach to elas tic-plastic fracture mechanics (EPFM) 
has been developed. Kardomateas and McClintock [42] indicate that 
ductile fracture occurs by hole growth, which depends on the combined 
history of strain, stress, and rotation near the crack tip. Fully 
plastic fracture mechanics (FPFM) can thus not be described by a single 
measure of failure. 

Several promising approaches to EPFM are being developed, however, 
which offer some insights into the process. These approaches generally 
consist of either a stress field analysis or an approximate global 
energy analysis. A stress field approach examines the stress field 
around the crack tip, as the K factor does in LEFM. An energy approach 
examines the energy release associated with crack initiation. Both 
approaches, unfortunately, have their limitations. 

Although it has not been accepted as a universal fracture 
criterion, the J- integral provides a means to determine an energy 
release rate for both LEFM and cases where plasticity effects are not 
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negligible [40], The J-integral is defined along a contour around the 
crack tip as the change in potential energy for a virtual crack 
extension. The U.S. Navy has adopted a J-based approach to fracture 
toughness characterizations of ductile materials [43], In general, 
though, the J-integral is unable to handle large amounts of stable crack 
growth accompanied by full plasticity because unloading occurs along a 
different s tress - strain curve than during loading, and also because of 
differing degrees of triaxiality with different geometries and low 
hardening. McClintock [10] further emphasizes that an energy approach 
to plastic crack propagation is global in nature and can not account for 
the different fracture process zone stress-strain distributions and the 
local mechanisms of crack growth. Crack initiation, however, is an 
especially local phenomenon and can be adequately described by an energy 
approach. 

Atkins [44] has developed a rigid-plastic fracture analysis which 
reasonably accounts for and simplifies the complex nature of fully 
plastic fracture (FPFM) in certain situations. In this approach, the 
relative magnitude of flow and fracture components is an important 
factor in determining the transition from plastic flow preceeding 
fracture to crack initiation and propagation. Although this approach is 
limited by its global nature, it offers a useful engineering analysis of 
FPFM, and is presented because of its promise. 

A limit load analysis offers a new approach to the problem of fully 
plastic crack propagation. By selecting realistic displacement fields 
it may be possible to define what loads cause deformation which could 
tend to promote crack propagation. 
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6 . 4 Rigid-Plastic Fracture Analysis 



A rigid-plastic simplification of EPFM may be applied to a variety 
of sheet tearing problems. Plastic fracture involves remote flow, which 
Atkins [44] incorporated into an expanded quasi-static elastic fracture 
energy equation in terms of the elastic strain energy A, the fracture 
toughness or specific (essential) work of fracture R, the total plastic 
work W p , the residual strain energy A , and kinetic energy H as: 

Pdu«dA + RdA + dW p + dA + dH (6.1) 

r 

where P is load, u is load-point displacement, and A is the crack area. 
In the specific case of dynamic tearing of ductile metal sheets the 
limiting crack speeds consistent with a quasi-static approximation are 
usually some 1 to 10 m/s [44] . The kinetic energy term may be equated 
to zero for most materials and cracked geometries. Cracks, whether 
surrounded by elastic or plastic fields, tend to propagate with their 
own crack opening angles. The rigid-plastic fracture equation can be 
further simplified in the case of extensive plastic flow, where it may 
be permissible to neglect dA, as: 

P du = R dA + dW p . (6.2) 
This equation has found success in Cotterell and Mai's [45] studies of 
fracture in ductile sheets, and in Atkin's studies of metal forming 
fracture such as cropping and can opening. These studies demonstrated 
that geometry- independent R values do characterise the fracture behavior 
and that plastic flow and fracture terms may be separated in the 
rigid-plastic fracture equation. Atkins cautions, however, that 
fracture toughness may be path dependent according to a McClintock- type 
analysis . 
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Atkins and others have examined crack propagation due to various 
loading modes. Experiments by Mai and Cotterell [46] in which the 
bending and unbending plastic work could be altered independently of the 
fracture work gave roughly the same values for fracture toughness. The 
experiment used two metals (low alloy steel Lyten and a 5251 alluminum 
alloy) with different work hardening indices in both two and three leg 
"trouser” tests to alter the plastic work independently of the fracture 
work. This demonstrates that the flow and fracture work components can 
sometimes be uncoupled and that the toughness so determined is a 
material property, and indicates that the rigid-plastic equation is 
valid for certain applications. 

From a global energy approach, crack propagation will occur if the 
combined flow and fracture associated with propagation requires a 
smaller work increment than flow alone. An increment of plastic flow 
may be written: 

dW p - d(r V) - r dV + V dr (6.3) 

where T is the total plastic work per unit volume and V is the volume 
undergoing plastic deformation. The plastic flow equation governs the 
external work required prior to fracture and applies with incremental 
changes dV and dW p since there is new material entering the fracture 
process zone as crack propagation occurs. By considering 
increments of external work, and the additional work required to 
propagate cracks, Atkins was able to describe stable crack 
propagation as : 




where n is the number of propagating splits [44]. The stability of 



77 



crack propagation is dependent upon boundary conditions and and loading 
as limited by the structure's comliance. 

This analysis was applied by Atkins [47] to the case of axial 
splitting of ductile metal tubes. The ends of the tubes were flared by 
being forced down onto conical mandrels, as shown in Figure 6.2. The 
number of starter cuts produced extremely interesting axial splitting 
behavior. With a small number of starter cuts, all cuts would propagate 
and some would bifurcate. With a larger number of starter cuts, some 
cuts would never develop into cracks. In all cases, the final number of 
propagating cracks was always between 8 and 12 for 50mm diameter and 
1.6mm wall thickness mild steel and aluminum tubes. Crack propagation 
was a stable process with a reasonably constant radius and all the 
splits having equal lengths. Atkins applied geometry and experience to 
the incremental plastic flow analysis in order to derive the 
relationship : 

n = (8tt/ 4T" )(a / R) r e 0 (6.5) 

where a is the effective stress, r is the radius of the tube, and e is 

u 

the average hoop strain. Atkins' method predicts n = 12, which 
indicates that the incremental plastic flow analysis can be successfully 
applied to some complex problems. 

6 . 5 Plastic Limit Load Fracture Analysis 

Crack propagation is encouraged by transverse strain, and ensuing 
necking and hole nucleation, in front of the crack tip. These strains 
can be produced by certain displacement fields within a fully plastic 
structure. The limit load for a crack of given length, with necking 
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neglected, can be found by examining the simple case of vertical up and 
down cracks extending from a central circular hole in a flat plate. 

An upper bound to the limit load is given by a displacement field 
which satisfies any displacement boundary conditions and gives no change 
in volume anywhere. The upper bound P can be found from: 

u b 

6W P = P Su p - f Y Sl p SV (6.6) 
ub J 

where W p is the plastic work, u p is the displacement under P in its 
direction, Y is the flow strength, € P is the equivalent plastic strain, 
and V is the volume. Experience with the crushing of initially deformed 
structures suggested separating the displacement field into two zones, 
as shown in Figure 6.3. The Zone II displacement field reflects the 
dominant shear strain evident during testing which places the crack 
under tension. Plastic work is equal to the applied load multiplied by 
the displacement caused by plastic deformation of the structure. the 
total displacement of the edge experiencing uniform compressive loading 
in the direction of the applied load is composed of in-plane compression 
and drop due to deflection. 

Plastic work is a function of stress and strain throughout the 
structure. For non-hardening plasticity, the plastic work occurs at the 
yield strength, which means that the flow stength Y can be considered 
equal to Since it is assumed that the thickness remains essentially 
constant during deformation, Eq . 6.6 can be rewritten as: 

W p = h JJ 8e P dxdy . (6.7 ) 

The integral must be evaluated separately for both zones and then summed 
to find the upper bound to the limit load. The equations are simplified 
by defining a displacement ratio of the maximum vertical deflection to 
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the maximum horizontal deflection. Since the equivalent plastic strains 
for this case are so complex, the plastic work is solved using a 
numerical integration technique. 

The results indicate that as the initial flare (dishing) angle 
increases the least upper bound decreases . Curves of normalized load 
versus displacement ratio for different initial angles of flare are 
presented in Figure 6.4. The normalized least upper bound decreases 
from 0.58 to 0.35 as the initial angle of flare increases from 0 to 45 
degrees. The displacement ratio associated with the normalized least 
upper bound also decreases with increasing initial angle of flare. 

Hodge [11] and Gordon and McCrum [12] performed uniform 
displacement limit load analyses of a flat square plate with a central 
circular hole of varying size, without any cracks. With a cutout ratio 
(c/b) of 40 percent, Hodge predicts an upper bound of 0.53 and a lower 
bound of 0.41, and Gordon and McCrum predict an upper bound of 0.60 and 
a lower bound of 0.47. The author's analysis describes a least upper 
bound of 0.58 for such a case. A structure will only fail in the 
assumed mode involving crack propagation if it does not first fail at a 
lower applied load according to a different mechanism. 

Although the author's approach is analytically correct, it does not 
provide any indication of the effect of an initial crack. It is not 
reasonable that structures will be unaffected by the presence of 
cracks, regardless of material properties. It seems likely that the 
least upper bound is too "loose,” and that the displacement field needs 
to be improved. Further upper bound limit load analysis of fractures 
could be conducted using a displacement field which specifically 
incorporates thinning in front of the crack tip, different u and v 
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Crack 



displacement fields, and various hole and crack sizes, 
propagation is promoted by a displacement field which produces a 
combination of strain and triaxial stress in front of the crack. A 
fracture criterion could then be assumed in the form of crack growth per 
unit thinning or displacement across the crack tip. 

6 . 6 Crack Propagation from a Circular Hole 

Crack propagation from a centrally located circular hole in a thin 
shell is a very complex problem. A J-based approach does not seem 
appropriate, Atkins' incremental plastic flow analysis can not be 
readily applied, and a limit load analysis does not provide significant 
illumination of the problem. Further research is required in the field 
of fully plastic crack propagation. 



81 



Figure 6 . 



Figure 6. 




: The three modes of loading which lead to crack growth [40] 




Axial splitting of ductile metal tube [47] 
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Figure 6.3: Quadrant of damaged plate and assumed displacement fields 
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7. EXPERIMENTS AND DISCUSSION 



7 . 1 Introduction 

An internal airblast produces a variety of damage modes examined in 
this paper. The combined structural effects of damage, which consist 
primarily of holing, dishing and tearing, have not been previously 
studied. A small program of experimentation, consisting of crushing 
seven different specimens, was therefore developed to gain qualitative 
insights into the interaction of various damage modes. 

7 . 2 Cumulative Structural Effect of Damage 

By considering the individual damage effects as a whole it is 
possible to paint a picture of the residual strength of structural 
plating within a ship's compartment after an internal explosion. 
Plastic limit load analysis of a flat plate indicates that ultimate 
strength decreases in proportion to increasing hole size. 
Elastic-plastic buckling and post-buckling analysis indicates that 
ultimate strength is little affected by hole size until the holing 
reaches half the plate width, and is essentially equal to the peak load 
defined by the effective width concept. The dishing of a flat plate by 
an explosive loading is related to impulse and plate geometry and 
material properties. An extension of the concept of small initial 
imperfections to a dished plate suggests that the ultimate strength is 
reduced by roughly 10 percent from that of a flat plate. Finally, the 
extensional collapse mode analysis indicates that the mean crushing load 
of shipboard structural plating is roughly 0.1 of the tensile yield 
load. Although the mean crushing load decreases with increasing hole 
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size, it is suggested that it is relatively unaffected by the extent of 
initial dishing. Fully plastic crack propagation has been examined, but 
its limiting effect can not be precisely described. Warship design 
incorporates the use of ductile structural plating and crack arrestors 
to negate the effect of initial cracks upon residual strength. 

Disregarding the potential effect of crack propagation, the 
theoretical combined effect of damage is shown in Figure 7.1 for 
shipboard structural plating with a length- to - thickness ratio of 250. 
Ultimate strength is described by the elastic -plastic out-of-plane 
failure mode presented in Chapter 3. The plastic in-plane failure mode, 
presented in Chapter 2, does not limit the peak load. Large-scale 
dishing, on the order of 5 to 10 plate thicknesses from explosive 
impulse, only reduces the ultimate load of an unperforated plate by 10 
percent. Following peak load, the mean crushing load ranges from 0.06 
to 0.13 for an unperforated plate, as presented in Chapter 5. Both 
ultimate and mean crushing strength are little affected by small holes 
(less than half the width) . While the cumulative effect of damage can 
significantly reduce the residual strength of structural plating, the 
plating retains roughly 40 to 50 percent of its original strength. The 
residual strength, even with extensive damage, is improved by support 
members which continue to provide clamped edge boundary conditions to 
structural plating. The compliance of the ship tends to restrict 
extensive crushing and provides support to the plating to allow it to 
continue to support normal service loads. 

7 . 3 Experimental Method 

Although not an exhaustive study, the experimentation was designed 
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to examine the impacts of different failure modes of damaged structure, 
and validate the analysis summarized in Chapter 7.2. The theoretical 
analysis is concerned with the case of an initially flat square plate. 
The simplest way to perform a uniaxial compressive loading of the plate 
was to construct a symmetrical box with four square side plates and two 
stiffened end plates. This arrangement produced a "natural" boundary 
condition at the side plate edges which maintained their right angle 
relationship during initial yielding, but did not prevent out-of -plane 
deformations. The side plates were essentially clamped at top and bottom 
and simply supported at the unloaded edges. 

A baseline "undamaged" specimen was built of four undeflected side 
plates. Two basic specimen types were designed to reflect different 
possible modes of lateral deflection produced by the airblast, as shown 
in Figure 7.2. In order to model the different possible modes of 
dishing (lateral deflection) from the initial undamaged condition, edge 
pieces were cut slightly larger and bent to shape such that the base 
dimensions and height of each specimen maintained the initial 
dimensions. Type A specimens represent an "Accordion bulge" which has 
caused the edges of the plate to deflect along with the center. Type B 
specimens represent a "Basketball bulge" which has not deflected the 
edges of the plate. Due to geometrical constraints of the supplied 
metal, Type A specimens were built with a peak deflection- to- thickness 
ratio of only 3.9, while the Type B specimens were built with a peak 
deflection- to- thickness ratio of 12.9. Three specimens were built of 
each type, with a central circular hole on all faces. 

The specimens were built of 1/16 inch thick mild steel (ASTM A-515 
grade 70) plates welded together by the TIG method. The standard ASTM 
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characteristics for this steel are presented in Table 7.1. The U.S. 
Navy uses ASTM A-515 grade 70 steel in construction of Navy boiler drums 
and headers. During fracture mechanics tests performed by the Navy, 
this carbon steel underwent substantial non-linear deformation followed 
by cleavage. It is in fracture mode transition at room temperature 
[48]. The U.S. Navy has recently adopted HSLA-80 (ASTM A-710 grade 
A-class 1 is typical) steel for use as structural plating in the hull 
and superstructure of warships. The standard ASTM characteristics for 
this micro-alloyed, fine-grain, low-carbon steel (as rolled and 
precipitation heat treated) are also presented in Table 7.1. 



Table 7.1: Tensile properties of ASTM A-515 grade 70 mild steel 

and ASTM-710 grade A-class 1 low-carbon steel [49] 



A-515 mild 

Composition (%) 


steel 

C < 0.27 , Mg 0.79-1.30 
P < 0.035 , S < 0.04 
Si 0.13-0.45 


, Ultimate Tensile Strength 


70-90 ksi [485-620MPa] 


, Yield Strength 


38 ksi [260 MPa] 


Maximum Thickness 


8 in [200 mm] 


A-710 HSLA 

Composition (%) 


steel 

C < 0.07 , Mg 0.40-0.70 
P < 0.025 , S < 0.025 
Si < 0.4 , Ni 0. 7-1.0 


* Ultimate Tensile Strength 


90 ksi (655 MPa) 


, Yield Strength 


80-85 ksi (550-585 MPa) 


Maximum Thickness 


3/4 in (20 mm) 



Hardness testing (Rockwell 30T) of the as-rolled specimen material 

2 

yielded Knoop numbers ranging from 112 to 120 kg/mm . Empirical 
analysis of these values indicated that the specimen material is softer 
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than that specified by ASTM characteristics. Since accurate tests of 
the yield strength were not conducted, however, the tensile yield load 
for the specimens was calculated with the specified yield strength 
value of 38 ksi. 

The initial undeflected square plate had dimensions of 5.5 by 5.5 
inches (140 by 140mm), for a length- to - thickness ratio of 90. A 
warship's structural plating has a length- to - thickness ratio, neglecting 
stiffeners, on the order of 200 to 300. Economic and time limitations 
dictated the selection of A-515 mild steel for constructing the 
specimens. Although this is a good ductile steel to represent the 

performance of HY-80 and HSLA-80 used in U.S. Navy warships, the 
specimens were relatively thicker and softer than actual structural 
plating. 

A simple undamaged specimen was crushed to serve as a baseline for 
comparison purposes. The first set of tests consisted of a specimen of 
each type, labeled Al and Bl respectively, with the lateral deformation 
pattern unique to each. Specimens Al and Bl had a 2-1/4 inch (57mm) 
central circular hole on all side plates, which represented a cutout 
ratio (c/b) of 41 percent. The second set of tests was similar to the 
first, but with the addition of a rosette of 8 cracks, 1/2 inch (12.6 
mm) long radiating from each hole, and labelled A2 and B2 . The holes in 
specimens A2 and B2 were the same size as in the first set of specimens. 
Crack tips were cut with hack-saw blades which had their sides ground 
down to ensure good initial notches for crack growth. The third set of 
tests was designed to evaluate the effects of changing hole size and 
initial crack formation. Specimen A3 had the same hole size as previous 
specimens, but was modified with a set of short 1/4 inch (6mm) vertical 
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cracks radiating from each hole. Specimen B3 was modified with a larger 
central circular hole of 3-3/8 inches (86mm), which represented a cutout 
of 61 percent, but had no initial cracks. Table 7.2 presents a 
tabulation of the initial conditions of the various specimens. The six 
"damaged" specimens represented a variety of initial damage effects. 



Table 7.2: Initial structural condition for all specimens 



Specimen 


Baseline 

Un- 

damaged 


Type A (Accordion) 


Type B (Basketball) 


Al 


A2 


A3 


B1 


B2 


B3 


w 

Dishing [ h ] 


0 


3.9 


3.9 


3.9 


12.9 


12.9 


12.9 


Cutout C b 3 


0 


0.41 


0.41 


0.41 


0.41 


0.41 


0.61 


Starter Cracks 
per Plate 


0 


0 


8 


2 


0 


8 


0 



The specimens were crushed on a Baldwin-Tate - Emery 60,000 psi 
machine to determine load versus displacement characteristics. The 
testing apparatus is shown in Figure 7.3. Uniform displacement was 
regulated by manual control of hydraulic flow. Load versus displacement 
data beyond ultimate load was only recorded by first unloading and then 
reloading slowly until the specimen resumed deflection. The testing 
machine had a nominal compliance of 1 inch per 200,000 lb. which was 
considered negligible for this experimentation. Loads were normalized 
with respect to the tensile yield load as described by Eq . 2.5: 

P *= a A (7.1) 

Y Y 

and displacements were normalized with respect to the original specimen 
height. Specimens were crushed to 13 to 50 percent of their original 
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length . 



7 . 4 Experimental Results 

All specimens exhibited generally similar elastic load- deflect ion 
curves, although ultimate strength and post-buckling behavior varied 
considerably. The original load-displacement data for all specimen 
crush tests is recorded in Apppendix A. All of the specimens 
demonstrated linear elastic behavior during the initial loading. Loss 
of strength was generally rapid after the peak load was reached. With 
continued crushing, however, the strength of each specimen leveled off 
and then increased as local folding zones collapsed. 

The baseline undamaged specimen supported a peak load of 20,600 lb. 
prior to collapsing, for a normalized ultimate load of 0.393. 
Normalized load versus displacement curves for all specimens are 
presented in Figures 7.4 and 7.5. The loads were normalized with 
respect to the tensile yield load , and the deflections were 
normalized with respect to the original baseline dimensions. The first 
holed specimens, Al and B1 , supported normalized ultimate loads of 0.247 
and 0.275 respectively. The buckling mode which developed in specimen 
B1 after the peak load was reached is shown in Figure 7.6. The second 
set of holed specimens with a rosette of cracks, A2 and B2 , supported 
normalized ultimate loads of 0.223 and 0.280 respectively. Only 
specimen A2 experienced crack propagation. The crack growth occurred 
only after the crack tips had been re-notched to eliminate the bluntness 
which had developed during the initial loading and after extensive 
crushing beyond the peak load. The buckling mode which developed in 
Specimen A2 after the peak load was reached is shown in Figure 7.7. The 
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limited crack propagation which developed was primarily Mode III tearing 
as shown in Figure 7.8. Specimen A3, with short vertical cracks, only 
supported a normalized ultimate load of 0.253. Although pinching 
(necking) was observed at the crack tips, specimen A3 experienced no 
crack propagation. Specimen B3, with larger circular holes, supported a 
normalized ultimate load of 0.248. The only visible weld failure 
occurred in specimen B3 during crushing beyond 40 percent of initial 
length. Table 7.3 presents a tabulation of the experimental results. 



Table 7.3: Ultimate and mean crushing load 

experimental results for all specimens 



Specimen 


Baseline 

Un- 

damaged 


Type A (Accordion) 


Type B (Basketball) 


A1 


A2 


A3 


B1 


B2 


B3 


Ultimate Load 
(lb.) 


20060 


12600 


11400 


12920 


14060 


14320 


12680 


Normalized 
Ultimate Load 


0.393 


0.247 


0.223 


0.253 


0.275 


0.280 


0.248 


Mean Crushing 
Load (lb.) 


8700 


9700 


9200 


9700 


5600 


5600 


5600 


Normalized 
Mean Crush Ld 


0.17 


0.19 


0.18 


0.19 


0.11 


0.11 


0.11 



7 . 5 Discussion of Results 

The experimentation gives some insight into the loss of strength 
associated with damage. The cumulative analysis of damage effect 
failure modes is applied to the specimen side plates and summarized in 
Figure 7.9. The experimental results are also plotted in Figure 7.9 for 
comparison purposes. Out-of -plane buckling caused all of the specimens 
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to collapse at lower loads than the tensile yield load. The effect of 
damage can, perhaps, be best seen by comparing damaged specimen strength 
characteristics with the baseline undamaged specimen. Although the 
experimentation consisted of just seven specimens, some trends or 
observations can be made about each of the different damage effects. 

Boundary conditions affect the plate buckling mode and strength 
characteristics. Figure 7.10 shows how the plate buckling coefficient 
K, which is directly related to the plate critical load, is affected by 
boundary conditions. The actual boundary conditions imposed upon the 
specimen side plates is illustrated by the clamped loaded edges - C 
condition shown in Figure 7.10. The plastic limit load and 
elastic -plastic analyses presented in Chapters 2 and 3 only examined 
simply supported edge boundary conditions. The critical buckling curve, 
shown in Figures 5.2 and 5.3, indicates that buckling of slender square 
plates of baseline specimen dimensions occurs at roughly 40 percent of 
peak load. Although the actual boundary conditions of the specimen 
plates should have given the plates greater strength than predicted 
analytically, out-of -plane buckling and peak load for the baseline and 
damaged specimens occurred at lower loads. 

All experimental results for ultimate and mean crushing load lie 
below the analytic solutions , as shown in Figure 7.9. Since the 
analytic solution for ultimate strength of an unperforated flat plate is 
known to be reliable, it seems likely that the low normalized ultimate 
strength for the baseline specimen was produced by using too high a 
yield strength value to determine the tensile yield load. Hardness 
testing indicated that the specimen material was softer than specified. 
A lower yield strength of 30 ksi brings the normalized experimental 
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results into line with the analytic solutions . Since the true yield 
strength is not known, however, it is useful to compare the results of 
crushing the damaged specimens with the baseline specimen. 

Increased holing clearly reduces residual strength. The behavior 
of specimens B1 and B3 during crushing describes the effect of increased 
holing. As the cutout ratio increased from 41 to 61 percent, with other 
initial damage effects constant, the normalized ultimate load decreased 
by roughly seven percent as compared to the baseline. These 
experimental results support the relative loss of strength with 
increased holing described by the elastic-plastic analysis. This 
indicates that although holing clearly reduces strength, it does not 
cause a rapid loss of strength . 

Dishing produces a significant loss of strength. It clearly biases 
the buckling mode towards its own shape. The experimentation ignores the 
effect of actual thinning, which would be associated with dishing, upon 
strength. While the undamaged baseline specimen buckled asymmetrically 
(opposite sides folded in or out together) , all damaged specimens 
buckled symmetrically with outward folds near mid-plate. Although the 
Type B specimens retained greater ultimate strength than the Type A 
specimens, their post-buckling mean crushing load was less. It seems 
reasonable that the Type B specimens retained greater ultimate strength 
because their edges were initially straight, while the Type A specimens 
were initially flared out. While the Type A specimens experienced 
horizontal folding along the mid-plane as predicted, the initial 
deflection (dished) shape of Type B specimens caused the folding zones 
to angle out from the central holes away from the mid-plane. In the 
post-buckling range, dished shape, holing, and slenderness may become 
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the dominant defining mean crushing load strength. 

A tentative quantitative estimate on the effects of dishing alone 
can be made by extrapolation of the experimental results. Assuming a 
linear reduction of strength with increased hole size, the peak load 
results from crushing of specimens B1 and B3 can be extrapolated to 
determine the loss of strength caused by dishing with no hole. A 
specimen with extensive initial dishing would thus have suffered a 16 
percent loss of strength from the baseline undamaged specimen, which is 
a slightly greater loss than predicted analytically. According to the 
elastic-plastic analysis, however, the results from crushing of 
specimens Al and B1 essentially represent the loss of strength due to 
dishing alone. A specimen with extensive initial dishing, but no 
holing, would thus have suffered a .30 to 32 percent loss of strength 
from the baseline undamaged specimen. The effect of actual thinning of 
structural plating during dishing caused by an explosion is considered 
negligible. There is unfortunately no way to determine the validity of 
this approach, or to apply it to other dished mode shapes, from the 
limited experimental data. 

Both initial hole size and deflection shape affect the mean 
crushing load. By assuming that the flow strength for A-515 mild steel 
equals 63 ksi, from Eq . 5.5 with o = 70 ksi, it is possible to apply 

u t s 

the extensional collapse mode analysis to the test specimens . The upper 
and lower bounds from the analysis describe a normalized mean crushing 
load ranging from 0.31 to 0.20, and from 0.24 to 0.17 for a specimen 
with a cutout ratio of 41 percent. The baseline undamaged specimen had 
a normalized mean crushing load of roughly 0.17 during crushing, which 
again suggests that the specimen material was softer than specified. 
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The extensional collapse mode analysis predicts a folding zone width of 
1.03 inches for the baseline specimen, from Eq. 5.7, which compares 
favorably to the observed width. Although the baseline specimen 
experienced repeated asymmetric folding, the damaged specimens 
experienced only one major folding pattern associated with the initial 
dished shape. Type A and B specimens, with a 41 percent cutout ratio, 
experienced different mean crushing loads of roughly 0.19 and 0.10 
respectively. Significant holing has the effect of allowing 
out-of-plane motion of the plating around the hole and thus reducing the 
amount of energy absorbed by extensional deformation of plate sections. 
The extensional collapse mode analysis does not describe a significant 
loss of mean crushing load with increased hole size until the cutout 
ratio exceeds roughly 60 percent, which is supported by the experimental 
results (specimen B3 versus Bl) . The initial angled folding zones 
associated with Type B specimen buckling reduced the range of folding 
(i.e. less than the complete folding presumed by the extensional 
collapse mode analysis) . It is suggested that holing and dishing 
combined can reduce the mean crushing load by roughly one- third as the 
dissipation of work becomes essentially quasi - inextensional deformation 
of hinge lines . 

Fracture is a potentially serious problem in any damaged 
structure, but was not a significant determining factor in the residual 
strength of the damaged specimens. The initial starter cuts seemed to 
have little effect upon strength, possibly because the fine cuts did not 
have any necking in front of the crack- tip as would have been present in 
natural fully plastic fracture. The presence of starter cracks had no 
consistent effect upon ultimate strength. The limited crack propagation 
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in specimen A2 was caused by tearing associated with complete folding 
and subsequent contact of buckles produced by extensive crushing. The 
mean crushing strength of the specimens was little affected by the 
presence of starter cuts. The thinning in front of the crack-tips 
indicate that transverse tensile strains were present, but insufficient 
to promote crack propagation. 

A ship's structural plating is relatively thinner and harder than 
the specimen side plates, and may experience quasi-static fracture more 
easily than the specimens did during testing. Fully plastic fracture in 
complex structures can not be described by current analytic techniques. 
Warship design, however, incorporates the use ductile materials to 
inhibit fracture, and stringers to serve as crack arrestors which limit 
any crack propagation. 

It is difficult to determine the comparative effects of damage 
between holing, dishing and cracking. The range of experimentation only 
covered two different damaged mode shapes, two hole sizes, and two crack 
patterns, as described in Table 7.2. Although further experimentation 
is required to quantify the individual and cumulative damage effects, it 
is possible to apply the lessons learned from this research. 

7 . 6 Residual Strength of a Warship After an Internal Explosion 

The DDG-51 Arleigh Burke Class destroyer is one of the U.S. Navy's 
most modern warships. The most demanding loading any ship 
generally experiences is located at its mid-section due to hogging and 
sagging conditions produced by the ocean. The DDG-51 is designed for a 
maximum hogging moment of 176,953 ft- tons and a maximum sagging moment 
of 87,883 ft-tons [50]. The maximum compressive loading condition of 
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2 

4.2 tons/in on the midships deck plating occurs during sagging. 

From a structural standpoint, the worst single antiship missile hit 
would be one which produces an internal explosion amidships. An 
internal airblast between frames 220 and 254, from a small antiship 
missile such as the EXOCET, could easily produce a hole with a cutout 
ratio of 50 percent and dishing on the order of 10 plate thicknesses. 
In terms of longitudinal strength, the deck plating above the damaged 
compartment can roughly be considered a 10 by 10 foot square plate of 
3/8 inch thick low-carbon steel. The worst condition sagging produces a 
load of 378,000 lb. on the plating, for a normalized load of 0.105. 
From the effective width concept, the undamaged structural plating can 
support a normalized ultimate load of 0.115, neglecting the positive 
effect of support members. Based upon the theoretical and experimental 
analysis presented in this paper, the ultimate strength of the damaged 
structure would be decreased by roughly 25 percent due to dishing and 5 
percent due to holing, for a normalized ultimate load of roughly 0.08. 
The normalized mean crushing load, incorporating the effect of holing 
but neglecting dishing, would range from roughly 0.07 to 0.04. These 
strength characteristics indicate that the plating can not support 
extreme sagging conditions prior to ultimate load on its own, but could 
probably continue to support normal service loads before and after peak 
load. A local failure of the structural plating in one or more 
compartments does not mean that the ship will experience global failure. 
The ship's support members and structural compliance would offer 
continued support to the plating to prevent local collapse. 
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Figure 7.1: Summary of damage effect failure modes for shipboard 

structural plating with a length- to- thickness ratio of 250 






Figure 7.2: Specimen types with different simplified dishing (lateral 

deflection) modes 
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NORMALIZED LOAD 




Figure 7.3: Photograph of testing apparatus 




Figure 7.4: Graph of normalized load versus deflection for baseline and 

Type A (Accordion bulge) specimens 
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NOnUAUZCD LOAD 




Figure 7.5: Graph of normalized load versus deflection for baseline and 

Type B (Basketball bulge) specimens 




Figure 7.6: Photograph of specimen B1 shortly after ultimate load was 

reached 



100 





Figure 7.7: Photograph of specimen A2 shortly after ultimate load was 

reached 
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normalized load 




Figure 7.9: Summary of damage effect failure modes and experimental 

results for Gspecimen side plates with a length- to - thickness ratio of 90 




Figure 7.10: Graph of plate buckling coefficient versus aspect ratio 

for different boundary conditions [18] 
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8. CONCLUSION 



Shipboard structural plating can continue to support normal service 
loading even after extensive damage. An internal explosion produces 
tearing, holing, and dishing of structural plating. The shipboard 
structure and damage effects were simplified in this paper in order to 
analyze the individual damage effects and determine how they interact. 
A small program of experimentation roughly validated the cumulative 
analysis of damage effect failure modes. 

Structural collapse is a function of mechanical properties, 
structural geometry, and collapse mode. (1) The ultimate strength of 
structural plating is strongly affected by slenderness. (2) The 
elastic-plastic analysis accurately describes the ultimate load for 
slender plates, and indicates that small holes (cutout ratio less than 
50 percent) have little affect. (3) Dishing on the order of 5 to 15 
plate thicknesses reduces ultimate strength by roughly 15 to 35 percent. 
(4) The buckling mode involves folding as decribed by the extensional 
collapse mode analysis. (5) The normalized mean crushing load, which 
describes the average compressive strength maintained after the ultimate 
load has been reached, is roughly 0.1 for undamaged shipboard structural 
plating, and ranges from 0.05 to 0.1 for significantly damaged plating. 
(6) Although crack propagation was limited during crushing of the 
damaged specimens, fully plastic fracture could be more significant in 
shipboard structural plating. Continued reserach is required to further 
quantify these conclusions and apply the results to improved ship 
designs . 
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8 . 1 Ideas for Continued Research 



This paper has presented the first qualitative theoretical and 
empirical analysis of the cumulative effect of damage upon structural 
plating. There are many avenues of continued research required to 
develop a comprehensive quantitative understanding of damage effects. 
Ideas for further research include: (1) determine the effect of ship's 
compliance upon plate boundary conditions, (2) determine the effect of 
stiffeners upon various damage modes, (3) investigate rectangular plate 
geometries, (4) determine the deflection of supporting members 
associated with plate geometries, (5) examine dynamic effects, (6) 
examine fully plastic fracture, (7) conduct extensive critical 
experimentation to determine quantitative relationship between damage 
effects. It is important that any testing properly model shipboard 
structural geometry (such as b/h) and material characteristics (such as 
^y/E) . The author suggests that determining the effcts of fully plastic 
fracture and shipboard structural compliance are particularly important. 
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APPENDIX A 

Load Versus Displacement Experimental Data 
from Crushing of Test Specimens 
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PHOTO 


—r c 


1 4720 




84. 75 


14010 




1 02 


1 3770 




1 15 


1 3220 
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10240 




199 


8720 




277 


5980 




712.5 


5150 




— c; -T ~r tr 


4620 




4 l 1 


4 0 3 0 


PHOTO 


570 


74 20 




oL5 . 5 


3250 




/ x 


3050 




6lu 


2970 




DEFLECT JON 


LOAD 




I rj * : _ *. ».< ) 


(LBS) 
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14 MAR 
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B3 - 



DEFLECTION LOAD 

( IN* 1000) (LBS) 



0 


(j 


7 


20 


12 


20 


21 


20 


31 


240 


35 


1 080 


38.5 


2170 


43.5 


3860 


47 


5280 


48. 75 


5950 


cr ^ cr 

vJ - . 4— 


6600 


cr —i 


7260 


54 


8080 


55.5 


8530 


57 


9080 


58 


9500 


60 


10090 


62. 5 


10640 


66.5 


1 1390 


6 . 5 


1 1 560 


71 


1 1 960 


74 


121 90 


S3 


1 255<J 


Q4 


12680 


1 06 


1 2 6 80 


1 26 


12440 


175 


1 0640 


1 97 


9880 


213 


9300 


2 2 6 


8580 


— C — > 
_ -J _ 


77 w 0 


278 


6250 


704 


5780 


338 


4 ”00 


44 1 


3 o B 0 


4 5° 
3 i 1 
550 


7-600 


718 0 


587 


7130 


C' — ~ 


3050 


-r- 4 


3000 


6 6 ” 


2970 


712 


2860 


^72 


2*9. d o 


756 


1850 


05* 


2° 10 


822 




850 


2 7 80 


887 


2 ”70 


943 


27 40 


96 6 


2740 



DEFLECTION 


LOAD 


t IN* 1 000) 


(LBS) 


1 006 


2780 


1032 


2790 


1088 


2790 


1 1 39 


2800 


1 1 86 


2850 


1259 


2850 


1 302 


29 1 0 


1 335 


2960 


1 380 


30 1 0 


1 420 


3080 


1478 


3 1 30 


1543 


3150 


1591 


3 1 70 


1 646 


3220 


1 706 


3230 


1732 


3280 


1 777 


3380 


1829 


3520 


1851 


3560 


1 Q 0 1 


3690 


1 942 


3720 


1984 


3880 


2030 


3980 


2077 


4060 
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